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1.  Introduction 


Breast-conserving  surgery  (BCS)  is  standard  treatment  for  breast  cancer.  However,  to  be  effective, 
surgery  must  remove  all  malignant  tissues.  Currently,  a  significant  fraction  of  BCS  patients  require 
additional  re-excision  surgery  due  to  positive  or  close  margins.  This  grant  aims  to  develop  a  system 
to  assess  tumor  margins  during  surgery,  with  the  eventual  goal  of  reducing  re-excision  surgery.  This 
new  technology  should  provide  accurate  information  about  the  status  of  the  margins  in  real-time,  and 
therefore  guide  decision-making  during  surgery  (Fig.  1).  These  technologies  are  designed  to  provide 
a  rapid  answer  such  that  the  surgeon  can  decide  on  the  spot  whether  to  remove  additional  tissue  to 
ensure  clear  margins. 


Figure  1.  Artistic  rendering  of  the  proposed 
approach.  A  Cerenkov  luminescence  imaging 
system  is  proposed  for  imaging  the  tumor 
margins  and  guiding  surgical  resection. 


The  impact  of  the  proposed  research  is  in  the 
avoidance  of  unnecessary  re-excision  breast 
surgeries,  since  it  is  hypothesized  that  the  proposed 
approach  will  help  visualize  compromised  margins 
directly  in  the  operating  room.  The  funded  project 
investigates  the  use  of  Cerenkov  luminescence 
imaging  for  providing  sensitive  and  specific 
information  about  the  cancer  during  surgery. 
Cerenkov  luminescence  is  the  natural  emission  of 
visible  light  from  radiotracers  such  as  FI  8- 
fluorodeoxyglucose  (FDG),  a  positron  emission 
tomography  (PET)  radiotracer  used  to  identify  highly 
metabolic  tumors  [1-3].  Since  FDG  is  retained  with 
high  specificity  by  breast  cancer  tumors,  it  is  hoped 
that  Cerenkov  luminescence  imaging  will  provide  a 
simple,  sensitive,  and  inexpensive  way  to  assess 
focal  accumulations  of  FDG  in  the  operating  room, 
thereby  allowing  highly  malignant  tissues  to  be 
identified  and  resected. 


2.  Body 

In  the  three  years  of  this  project,  the  following  milestones  have  been  achieved  and  the  following  tasks 
have  been  completed: 


Task  1.  Set-up  Cerenkov  luminescence  characterization  and  quantitation  platform. 

This  first  task  involved  the  development  of  an  imaging  system  for  visualizing  and  measuring  minute 
amounts  of  radioactive  molecules.  Two  systems  were  set  up  and  investigated.  The  first  system  was 
designed  for  imaging  very  small  samples  (e.g.  tissue  sections)  and  set  up  by  retrofitting  an  existing 
bioluminescence  microscope  (LV200,  Olympus)  available  to  us  (Fig.  2a).  The  LV200  is  designed  to 
image  bioluminescence  in  living  cells  rather  than  radioactive  tracers  such  as  FDG.  Hence,  we 
modified  some  of  the  hardware  and  software  on  that  system  to  allow  for  imaging  of  radioactive 
tracers.  The  LV200  is  an  inverted  microscope,  comprising  an  electron-multiplying  charge-coupled 
device  (EM-CCD)  camera,  a  high-numerical-aperture  lens  (40X,  1.3  NA),  and  live  cell  imaging 
capabilities.  It  is  equipped  with  various  filters  for  fluorescence  imaging.  The  entire  microscope  is 
shielded  from  room  lights  by  a  light-tight  enclosure.  Its  hardware  specifications  allow  it  to  resolve  very 
low  amounts  of  light,  such  as  those  produced  by  Cerenkov  from  radiotracers. 
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Figure  2:  (a)  Sensitive  microscope  for  measuring  Cerenkov  luminescence  and  beta  radiation  at  the  cell  level,  (b) 
Imaging  system  for  measuring  Cerenkov  luminescence  in  larger  samples,  (c)  FDG  drop,  imaged  with  the 
microscope.  The  color  overlay  represents  the  Cerenkov  emission  while  the  grayscale  image  is  the  brightfield.  (d) 
Image  intensity  in  a  region  of  interest  for  a  drop  of  FDG  as  the  activity  undergoes  radioactive  decay,  (e)  Brightfield 
image  of  MDA-MB-231  breast  cancer  cell,  (f)  Distribution  of  FDG  in  the  same  cells  measured  using  a  novel 
technique,  involving  placing  a  scintillator  plate  in  close  contact  to  the  cells,  (g)  Overlay,  (h)  Cerenkov  signal 
(color)  emanating  from  a  syringe  filled  with  FDG,  acquired  with  a  custom  imaging  system,  (i)  Emission  spectrum 
of  Cerenkov  luminescence,  measured  using  optical  bandpass  filters. 

The  second  system  was  designed  for  imaging  larger  specimen,  such  as  primary  tumors  or  small 
animals.  The  system  was  custom-built  and  is  composed  of  a  large  light-tight  enclosure,  a  moveable 
stage,  a  set  of  LEDs  for  taking  photographs,  and  a  sensitive  EM-CCD  camera  capable  of  measuring 
Cerenkov  luminescence  (Fig.  2b).  Software  was  developed  for  operating  the  system,  including  the 
capability  to  move  the  stage,  activate  the  LEDs,  and  take  images  with  the  camera. 

A  spectrograph  attachment  was  also  set-up  for  acquiring  spectral  measurements  of  Cerenkov 
luminescence. 

Task  2.  Implement  combined  beta  /  Cerenkov  detection  capabilities 

As  a  preliminary  experiment,  a  small  drop  of  FDG  (2  pCi)  was  deposited  in  a  glass-bottom  dish  and 
imaged  with  the  modified  LV200  microscope.  Cerenkov  luminescence  from  FDG  was  clearly  visible, 
showing  that  the  instrument  has  enough  sensitivity  to  detect  Cerenkov  luminescence  (Fig.  2c).  As  the 
radioactivity  decayed,  further  images  were  taken  to  establish  the  lowest  amount  of  radioactivity 
detectable  with  this  set-up,  which  was  found  to  be  approximately  1  nCi  (Fig.  2d). 

A  question  that  arose  from  these  preliminary  measurements  was  what  is  the  smallest  resolvable 
object  with  this  approach.  More  specifically,  could  single  cells  be  resolvable  using  Cerenkov 
luminescence?  To  answer  this  question,  breast  cancer  cells  (MDA-MB-231)  were  cultivated  in  glass- 
bottom  dishes  and  incubated  with  FDG  (300  pCi).  However,  no  Cerenkov  signal  could  be  resolved. 
The  experiments  were  repeated  several  times.  A  computer  simulation,  performed  later,  showed  that 
indeed,  a  single  cell  could  not  be  resolved  with  Cerenkov  luminescence  because  (1)  the  beta  particle 
can  travel  up  to  100  pm  in  tissue,  spreading  the  signal  over  a  substantial  volume;  and  (2)  the  signal- 
to-noise  ratio  remains  low  because  of  the  small  number  of  photons  emitted  by  Cerenkov 
luminescence  (on  average,  3  photons  per  decay). 
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Because  of  these  limitations,  we  looked  for  an  alternative  approach  that  could  visualize  uptake  of 
FDG  in  single  cells.  This  requirement  is  critical  to  our  research  goal,  that  is,  ensuring  that  no  cancer  is 
left  after  surgery.  To  solve  the  two  aforementioned  issues,  we  investigated  the  use  of  an  indirect  beta 
detection  scheme.  We  hypothesized  that  an  inorganic  scintillator  could  (1)  reduce  the  range  of  the 
beta  particle  due  to  the  high  density  and  high  atomic  number,  and  (2)  increase  the  signal-to-noise 
ratio  by  amplifying  the  number  of  optical  photons  (on  average,  a  beta  particle  can  produce  10,000 
photons  in  suitable  scintillator  materials).  We  investigated  CdW04  as  a  suitable  material  for  detecting 
low  amounts  of  radiotracers  such  as  FDG  in  cells.  To  maintain  close  contact  between  the  cell  and  the 
scintillator  plate,  we  cultivated  breast  cancer  cells  directly  on  the  scintillator  plate.  After  incubation 
with  FDG,  the  cells  were  imaged  with  the  modified  microscope.  The  images  showed  focal 
accumulation  of  FDG  in  single  cells,  co-localized  with  the  cells  position  (Fig.  2e-g).  This  new 
technique,  which  we  call  radioluminescence  microscopy,  is  unique  and  a  significant  advance: 
currently,  no  other  approach  can  measure  and  quantify  the  uptake  of  a  radiotracer  in  single  living 
cells.  These  new  methods  and  results  were  published  in  PLOS  One  (Appendix  1). 

Task  3.  In  breast-tissue-  and  tumor-mimicking  phantoms,  characterize  light  and  beta  propagation  and 
detection 

To  further  characterize  Cerenkov  luminescence  in  biological  tissues,  we  developed  a  computer 
simulation  framework  to  model  the  propagation  of  Cerenkov  photons  in  turbid  media  (i.e.  with  high 
scattering  and  absorbing  properties).  The  simulation  framework  used  Monte  Carlo  methods  to 
simulate  the  history  of  individual  photons.  For  accurate  results,  billions  of  photons  need  to  be 
simulated,  which  can  require  very  long  run  time  on  large  computer  clusters.  Therefore,  we  use  two 
different  approaches  to  accelerate  computer  simulations  and  make  them  practical.  The  first  approach 
used  a  “roulette”  mechanism  to  keep  photons  alive  even  after  they  have  been  absorbed.  In  this 
scheme,  the  simulation  models  packets  of  photons  rather  than  individual  photons,  which  help  reduce 
the  computational  requirements  tremendously.  Furthermore,  we  implemented  a  distributed 
computation  software  package  to  run  these  simulations  in  a  cloud-computing  environment.  Cloud 
computing  is  a  new  paradigm  that  allows  users  of  computing  resources  to  rent  these  resources  on 
demand,  for  a  given  amount  of  time,  rather  than  purchase,  install,  and  maintain  large  computer 
clusters  on  site.  Because  we  did  not  have  access  to  large  computer  clusters,  cloud  computing 
allowed  us  to  run  large  optical  simulations  on  240  nodes  in  less  than  two  hours.  Such  simulations 
would  have  taken  more  than  a  week  on  a  single  computer.  These  new  methods  and  results  were 
published  as  a  research  article  in  the  Journal  of  Biomedical  Optics  (Appendix  2). 

Using  this  simulation  tool,  we  were  able  to  better  characterize  the  propagation  of  Cerenkov  light  in 
biological  tissues.  Using  optical  properties  specific  to  breast  adipose  tissue,  we  simulated  the 
emission  of  photons  from  a  point  source,  using  the  wide  spectrum  characteristic  of  Cerenkov 
radiation.  The  depth  of  the  emission  source  was  varied  from  0  to  1  cm.  The  photon  detector  was 
assumed  to  be  ideal.  We  found  that  Cerenkov  measurements  could  be  performed  up  to  a  depth  of 
approximately  5  mm,  which  is  suitable  for  assessing  disease  margins. 

The  optical  emission  spectrum  of  Cerenkov  luminescence  was  also  measured  using  2  ml  of  FDG  in  a 
syringe  using  our  custom-built  imaging  system  (Fig.  2h).  The  wavelength  was  found  to  be 
proportional  to  the  inverse  of  the  wavelength  squared,  which  confirms  the  theoretical  observations  of 
previous  investigators  (Fig.  2i). 


While  year  1  was  focused  on  modeling  of  light  diffusion  using  computer  models,  in  year  2  we 
performed  experiments  intended  to  demonstrate  Cerenkov  luminescence  imaging  in  breast-mimicking 
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phantoms.  We  investigated  several  materials  that  have  optical  properties  similar  to  breast  tissue,  and 
found  the  following  mixture  to  have  the  best  properties:  2%  agarose  for  mechanical  strength,  1% 
intralipd  for  optical  scatter,  17  uM  hemoglobin  for  optical  absorption,  and  15  mml  NaN3  for  antibiotic 
and  antifungal  properties.  Using  high-melt  agarose,  we  first  fabricated  a  cancerous  lesion  that 
contained  200  uCi  of  FDG.  The  mixture  was  poured  in  a  1  cm  diameter  spherical  mold,  which  we 
fabricated  as  well.  We  let  the  mixture  cool  down  and  solidify.  We  then  embedded  the  lesion  in  a 
hemispherical  breast-shaped  mold  filled  with  low  melt  agarose.  The  cancerous  lesion  is  therefore 
entirely  surrounded  by  “normal”  tissue.  A  color  photograph  of  the  phantom  is  shown  in  Fig.  3A. 


We  then  placed  the  phantom  in  a  custom  Cerenkov  imaging  system.  The  system  was  built  in  Year  1 
and  was  reported  in  the  2011  summary  report.  First,  we  took  a  Cerenkov  luminescence  image  of  the 
intact  phantom.  No  Cerenkov  signal  was  detectable  because  the  lesion  was  too  deep.  After  incision 
of  the  phantom,  we  were  able  to  visualize  strong  Cerenkov  signal  emanating  from  the  exposed  lesion 
(Fig.  3B).  We  used  this  information  to  surgically  remove  a  piece  of  phantom  tissue,  which  we 
pelieved  was  the  lesion.  By  taking  another  Cerenkov  image,  we  were  able  to  verify  that  the  excised 

sample  was  indeed  malignant  due  to  the  high  Cerenkov 
signal  (Fig.  3C).  We  also  found  that  the  lesion  resection 
had  not  been  complete:  there  was  significant 
luminescence  coming  from  the  tumor  bed.  We  therefore 
performed  a  second  resection,  after  which  no  further 
Cerenkov  luminescence  could  be  observed  in  the  tumor 
bed  (Fig.  3D).  We  then  cut  the  phantom  open  to  ensure 
that  the  entire  lesion  had  been  removed.  Those 
experiments  were  repeated  three  times  to  ensure 
statistical  significance. 

This  experiment  helped  show  that  Cerenkov 
luminescence  imaging  of  FDG  is  a  promising  direction 
for  guiding  breast  cancer  surgery.  However,  we 
encountered  a  few  difficulties.  According  to  Dr  Wapnir, 
an  experienced  breast  surgeon  who  provided  her 
expertise,  the  consistency  and  mechanical  strength  of 
the  agarose-based  phantoms  was  significantly  different 
from  human  breast  tissue.  Furthermore,  we  noticed  that 
some  of  the  FDG  was  able  to  diffuse  away  from  the 
primary  lesion  and  into  the  surrounding  matrix.  FDG  does  not  diffuse  in  vivo  because  it  is  retained 
within  the  cells.  We  will  investigate  other  materials  that  have  optical  and  mechanical  similar  to  breast 
tissue.  In  the  meantime,  we  will  also  report  the  experiments  described  here  in  a  biomedical  optics 
journal. 


cancer\  Margjns 
removed  dear 


Figure  3.  Simulation  of  breast  cancer 
surgery  guided  by  Cerenkov  luminescence 
imaging  using  a  phantom 


Task  4.  Design  algorithm  for  estimating  signal  depth  based  on  combined  spectral  Cerenkov  and  beta 
scintillation  imaging 


We  investigated  penetration  depth  for  Cerenkov  imaging  and  direct  beta  detection  with  a  scintillator. 
Overall,  due  to  the  positron  range,  direct  beta  detection  with  a  scintillator  was  limited  to  sources  of 
radiation  less  than  200  urn  deep.  This  makes  direct  beta  detection  most  suitable  for  probing  shallow 
tissues,  such  as  residual  cancer  in  the  surgical  field.  Cerenkov  luminescence  was  detected  up  to  a 
depth  of  5  mm  (in  tissue-mimicking  material,  given  100  uCi  of  activity). 


We  found  that  one  of  the  advantages  of  the  direct  detection  approach  was  that,  due  to  the  physics  of 
scintillation,  more  photons  were  produced  compared  to  Cerenkov  luminescence.  To  quantify  this 
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effect,  we  measured  the  luminescence  intensity  using  both  methods  (Fig.  4A)  and  found  30X 
stronger  luminescence  using  a  LSO  scintillator. 

Figure  4.  (A)  Luminescence 
intensity  shown  for  Cerenkov 
luminescence  and  LSO 
scintillator.  The  scintillator 
provides  ~30X  signal 
enhancement.  (B) 

Luminescence  enhancement 
shown  for  a  flexible 
scintillator  film,  applied 
directly  on  an  excised  tumor 
specimen. 

101  io°  io_1  io"2  icr3 

FDG  activity  (pCi) 

This  effect  suggests  that  we  can  use  a  scintillator  to  amplify  low  concentrations  of  FDG  in  malignant 
tissue.  To  this  purpose,  we  fabricated  thin,  flexible  alginate  films  that  embed  nano-sized  scintillator 
powders  (here,  CaF2:Eu).  To  evaluate  these  scintillating  films,  we  wrapped  them  around  mouse 
xenograft  tumors  that  were  excised  1h  after  FDG  injection.  We  found  that  the  film  lead  to  a  97% 
increase  in  luminescence  intensity  (Fig.  4B).  We  are  planning  to  use  these  films  to  image  FDG 
uptake  in  excised  samples  with  higher  sensitivity. 

Task  5.  Design  algorithm  for  tomographic  imaging  of  tracer  distribution  in  a  tissue  sample  from 
several  projection,  using  the  combined  spectral  Cerenkov  and  beta  scintillation  signal 
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In  the  past  two  years,  several  methods  for  performing  Cerenkov  tomographic  reconstructions  have 

been  presented  [4,5].  Since  these  algorithms 
are  well  adapted  to  the  final  goal  of  this 
project,  we  did  not  develop  new  reconstruction 
algorithms  for  tomographic  reconstruction.  We 
have  evaluated  one  of  these  algorithms  on 
phantom  datasets  and  have  achieved 
satisfying  reconstructions. 


Figure  5.  Multimodal  imaging  of  two  mice  bearing 
MDA-MB-231  xenograft  tumors.  (A)  PET/CT  imaging; 
(B)  Cerenkov  luminescence  imaging;  (C) 
Autoradiography;  and  (D)  Radioluminescence 
microscopy.  Arrows  show  the  two  tumors. 


radioluminescence  microscope)  developed  in  year 


excision  of  the  tumors  (Fig.  5B).  The 
specimens  were  frozen  and  sectioned,  then 
imaged  with  standard  autoradiography  (Fig. 
5C).  Last,  slices  from  the  same  specimen 
were  also  imaged  using  the  high-spatial- 
resolution  imaging  system  (called 
1  of  this  project.  A  full  report  on  this  imaging 


Task  6.  Determine  the  biodistribution  of  FDG 
in  breast-cancer  mouse  xenografts  at  various 
time  points,  and  correlate  with  Cerenkov 
imaging,  autoradiography  and  pathology 

In  a  preliminary  study,  we  implanted  two  mice 
with  MDA-MB-231  human  xenografts.  The 
mice  were  injected  with  FDG,  and  imaged 
using  standard  PET/CT  (Fig.  5A).  After 
euthanasia,  the  mice  were  imaged  with 
Cerenkov  luminescence  before  and  after 
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system  was  recently  published  [6]  and  is  attached  as  Appendix  1.  The  results  of  this  experiment 
suggest  that  Cerenkov  luminescence  is  best  adapted  for  surveying  the  surgical  field,  because  it 
allows  the  visualization  of  a  large  area  (e.g.  Fig.  1).  The  radioluminescence  microscope  is  more 
suitable  to  analyzing  small  excised  specimen  with  high  spatial  resolution.  It  can  simultaneously 
visualize  radionuclide  probe  distribution  and  tissue  morphology  for  accurate  assessment  of  tumor 
margins  with  very  high  spatial  resolution.  FDG  uptake  in  single  cancer  cells  can  be  visualized  to 
ensure  that  the  outer  rim  of  tissue  on  the  excised  specimen  (i.e.  the  margin)  is  cancer-free. 

In  year  2  of  this  project,  we  have  also  developed  a  new  image  acquisition  method  for  the 
radioluminescence  microscope  system  that  is  more  accurate.  Rather  than  using  a  long  exposure  (e.g. 
5  min),  the  new  method  takes  thousands  of  images,  each  with  a  short  exposure  (e.g.  100  ms).  This 
allows  single  molecules  of  FDG  to  be  visualized  as  they  decay.  This  particle  counting  scheme  has 
two  advantages:  (1)  the  spatial  resolution  is  improved  because  each  individual  radioactive  decay 
event  can  be  localized  with  greater  precision;  and  (2)  the  images  have  higher  quantitative  accuracy 
because  each  decay  event  is  counted  as  “one”  and  not  a  random  quantity.  This  new  imaging  scheme 
was  described  in  a  new  journal  article  published  in  the  Journal  of  Nuclear  Medicine  (Appendix  3). 

Task  7.  Develop  bedside  tissue  characterization  platform 

As  a  first  step  toward  translation,  we  developed  an  endoscopic  Cerenkov  imaging  system  that  is  more 
applicable  to  clinical  human  imaging  [7].  This  system  is  based  on  a  fiber  bundle  connected  to  a 
sensitive  image-intensified  CCD  camera  via  an  optical  coupler  (Fig.  6).  This  fiber  bundle  system  is 
flexible  and  can  be  used  to  visualize  Cerenkov  luminescence  signal  at  various  body  sites,  including 
inside  cavities  such  as  the  gastrointestinal  tract.  In  order  to  achieve  sensitivity,  the  ICCD  camera  is 
run  in  single  photon  mode  by  thresholding  individual  photon  packets.  We  were  the  first  group  to  report 
the  possibility  of  imaging  Cerenkov  luminescence  through  an  endoscopic  device  in  an  article 
published  in  the  Journal  of  Nuclear  Medicine  (Appendix  4). 


Figure  6.  Fiberscopic  prototype  system  for  clinical  imaging  of  tissue  samples.  The  system  (shown  on  the 
left)  comprises  a  flexible  endoscope  that  can  be  brought  close  to  the  sample,  the  system  was  sued  to  image 
a  glioma  xenograft  mouse  model  (top  right)  and  was  compared  to  a  commercial  small-animal  imaging 
system  (bottom  right). 

In  order  to  further  evaluate  the  performance  of  the  endoscopy  system  in  living  system,  we  developed 
murine  models  of  human  cancer.  We  selected  a  C6  glioma  cancer  model  for  this  study  because 
cancer  xenografts  based  on  this  model  are  more  infiltrative  and  more  representative  of  actual  breast 
cancer,  at  least  from  a  physical  perspective.  The  C6  glioma  model  is  more  likely  to  result  in  residual 
cancer  after  surgery,  which  makes  it  an  ideal  model  for  evaluating  Cerenkov  luminescence  imaging. 
Five  mice  bearing  G6  xenografts  were  injected  with  FDG  (1  mCi)  and  imaged  with  both  the 
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fiberscopic  system  and  a  commercial  optical  imaging  system  (IVIS  200).  The  mice  were  sacrificed 
prior  to  starting  surgery.  The  new  fiberscopic  probe  showed  that  it  was  possible  to  visualize  Cerenkov 
luminescence  emanating  from  FDG  in  the  tumor.  The  boundaries  of  the  xenograft  tumor  were  clearly 
delineated  and  displayed  some  heterogeneity.  The  image  quality  obtained  with  this  set-up  was 
comparable  to  that  achievable  with  a  commercial,  small  animal  imaging  system.  However,  it  should 
be  noted  that  small-animal  imaging  systems  do  not  offer  the  same  flexibility  a  the  new  fiberscopic 
system,  which  can  be  bent  to  access  internal  cavities.  In  conclusion,  we  were  able  to  design  and  build 
a  prototype  imaging  system  that  could  one  day  be  used  in  the  clinic  for  assessment  of  surgical 
margins  of  resection. 

Training:  As  part  of  the  comprehensive  training  plan,  I  have  attended  various  events  to  increase  my 
knowledge  of  breast  cancer  research.  While  conducting  the  research,  I  have  learned  how  to  conduct 
the  following  experimental  procedures:  cancer  tissue  culture,  live  cell  microscopy,  modeling  of  light  in 
biological  tissues,  use  of  radiopharmaceuticals  such  as  FDG,  optical  instrumentation,  EM-CCD 
cameras,  and  scintillator  physics.  Furthermore,  during  the  course  of  the  research,  I  was  able  to 
interact  with  various  researchers  with  a  very  diverse  set  of  background  and  skills.  The  interactions 
with  my  co-mentor  Dr.  Irene  Wapnir  (breast  cancer  surgeon)  were  particularly  useful  and  productive.  I 
was  given  the  unique  opportunity  to  observe  various  breast  cancer  procedures.  I  also  took  advantage 
of  the  exceptional  educational  resources  available  to  postdocs  on  the  Stanford  campus.  For  instance, 
I  took  several  classes  highly  relevant  to  this  proposal  research  aims,  including  “Biophotonics:  Light  in 
Medicine  and  Biology”  (EE  331),  the  Stanford  Comprehensive  Cancer  Research  Training  Program 
(see  attached  certificate),  and  “Multi-modality  molecular  imaging”  (BIOE  222).  Furthermore,  I 
attended  a  number  of  conferences  on  molecular  imaging,  medical  physics,  biophysics,  and 
personalized  medicine. 

3.  Key  Research  Accomplishments 

•  Development  and  characterization  of  two  platforms  for  imaging  Cerenkov  luminescence  both  in 
small  and  large  samples. 

•  Cerenkov  luminescence  was  observed  from  decay  of  FDG  in  millimeter-sized  samples,  but 
could  not  resolve  single  cells  due  to  poor  spatial  resolution  and  efficiency. 

•  A  new  approach,  called  radioluminescence  microscopy,  was  developed  to  measure  and 
quantify  radiotracer  concentration  in  single  cells.  The  new  approach  use  a  scintillator  plate  to 
improve  spatial  resolution  and  efficiency  compared  to  Cerenkov  luminescence. 

•  The  new  approach  was  used  to  investigate  uptake  of  FDG  by  single  breast  cancer  cells. 

•  A  new  computer  simulation  framework  was  developed  for  accurately  modeling  emission  and 
propagation  of  Cerenkov  luminescence  in  biological  tissues. 

•  Phantom  study  simulating  breast  cancer  surgery  guided  by  Cerenkov  luminescence  imaging 

•  Investigation  of  Cerenkov  imaging  at  different  depths 

•  New  method  developed  to  image  single  radionuclide  decay  events  with  cellular  resolution, 
resulting  higher  image  quality  and  accuracy  for  radioluminescence  microscopy 

•  First  demonstration  of  radioluminescence  microscopy  for  tissue  samples,  such  as  tumor 
specimen  excised  during  breast  cancer  surgery 

•  Imaging  of  breast  cancer  xenograft  models  using  multiple  modalities:  PET/CT,  Cerenkov 
luminescence  imaging,  autoradiography,  and  radioluminescence  microscopy 

4.  Reportable  Outcomes 

•  Four  journal  articles  published,  four  currently  under  preparation  /  review 
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Pratx  G  &  Xing  L,  "Monte  Carlo  simulation  of  photon  migration  in  a  cloud  computing 
environment  with  MapReduce",  J.  Biomed.  Opt.  16(12)  pp.  125003,  Dec.  2011 

Pratx  G,  Chen  K,  Sun  C,  Martin  L,  Carpenter  CM,  Olcott  PD  &  Xing  L,  "Radioluminescence 
microscopy:  Measuring  the  heterogeneous  uptake  of  radiotracers  in  single  living  cells",  PLOS 
One  7(10),  e46285  (2012) 

Liu  HG,  Carpenter  CM,  Jiang  H,  Pratx  G,  Sun  C,  Buchin  MP,  Gambhir  SS,  Xing  L,  &  Cheng  Z, 
"Intraoperative  imaging  of  tumors  using  Cerenkov  luminescence  endoscopy:  A  Feasibility 
Experimental  Study",  J.  Nucl.  Med.  53(10),  pp.  1579-1584  (2012) 

G.  Pratx,  K.  Chen,  C.  Sun,  C.  Carpenter,  &  L.  Xing,  “High-resolution  radioluminescence 
microscopy  of  FDG  uptake  by  reconstructing  the  beta  ionization  track”,  J.  Nucl.  Med.,  54(10) 
pp.  1841 -1846  (2013) 

CM  Carpenter,  C  Sun,  G  Pratx,  K  Chen,  SS  Gambhir,  L  Xing,  “Multiscale  PET/Cerenkov 
Image-Guidance  to  aid  Oncological  Resection  in  an  Invasive  Prostate  Tumor  Mouse  Model”, 
submitted  to  Radiology 

C.M.  Carpenter,  X.  Ma,  H.  Liu,  Conroy  Sun,  G.  Pratx,  J.  Wang,  S.S.  Gambhir,  L.  Xing,  Z. 
Cheng,  “Improved  Cerenkov  Molecular  Sensitivity  with  Beta  (minus)  Emitting  Radiotracers”, 
submitted  to  J.  Nucl.  Med. 

S.  Turkcan,  Andrew  Ceballos  &  G.  Pratx,  “First  demonstration  of  10  pm  resolution  imaging  of 
18F-FDG”,  under  preparation 

A.  Natarajan,  S.  Turkcan  &  G.  Pratx,  “Multiscale  and  multimodality  imaging  of  radionuclide 
probes”,  under  preparation 

•  Three  conference  presentations  (oral) 

Pratx  G,  “Pushing  PET  imaging  to  the  cellular  level:  Development  of  a  radioluminescence 
microscope”,  AAPM  annual  meeting,  August  2013 

Carpenter  CM,  Sun  C,  Pratx  G  &  Xing  L,  “Multiscale  PET/Cerenkov  image-guided  surgery: 
Demonstration  in  an  invasive  mouse  tumor  model”,  AAPM  Annual  Meeting,  2013 

Pratx  G,  Chen  K,  Sun  C,  Martin  LM,  Carpenter  CM  &  Xing  L,  "A  new  method  to  measure 
radiotracer  uptake  in  single  living  cells",  oral  presentation,  WMIC,  2012  (Selected  as  a  Young 
Investigator  Award  semi-finalist) 

•  One  US  patent  application 

Imaging  the  Heterogeneous  Uptake  of  Radiolabeled  Molecules  in  Single  Living  Cells,  US 
patent  application,  13/492,606 

•  One  travel  award  to  attend  WMIC  in  Dublin  (2012) 

•  One  class  taught  as  lead  instructor  (Mechatronics  system) 

•  Three  classes  taken  (Photonics  in  Biomedicine,  Molecular  Imaging,  and  comprehensive  cancer 
research  training  program) 

•  On  April  1st  2013,  started  tenure-track  position  in  Radiation  Oncology  at  Stanford  University 


11 


5.  Conclusions 


The  first  phase  of  this  project  showed  that  Cerenkov  luminescence  could  be  easily  detected  from 
FDG,  a  radiotracer  that  is  preferentially  taken  up  by  malignant  tumors.  We  were  not  able,  however,  to 
detect  this  effect  from  single  isolated  cells.  To  solve  this  issue,  we  developed  a  new  technique, 
termed  radioluminescence  microscopy,  that  can  visualize,  measure,  and  quantify  uptake  of 
radiotracers  in  single  cells.  Unlike  Cerenkov  imaging,  radioluminescence  microscopy  uses  a 
scintillator  plate  to  amplify  the  ionization  signals.  Currently,  much  research  effort  is  dedicated  to  the 
characterization  of  the  parameters  of  single  cells,  as  opposed  to  large  cell  populations.  This  new 
technique  will  fill  a  critical  need  and  reveal  the  mechanisms  that  govern  the  uptake  of  radiotracers  by 
single  cells.  Furthermore,  this  significant  advance  demonstrates  that  it  may  be  possible  to  detect 
individual  cancer  cells,  which  should  provide  the  most  sensitive  assessment  of  residual  disease. 

In  the  second  phase  of  this  project,  we  have  further  developed  a  fiberscopic  system  designed  to 
image  FDG  uptake  during  breast  cancer  surgery.  The  strategy  for  assessing  residual  disease  during 
breast  surgery  consists  of  three  steps:  before  surgery,  standard  PET/CT  is  used  to  evaluate  how 
extensive  the  disease  is,  based  on  the  distribution.  Then,  during  surgery,  Cerenkov  luminescence  is 
used  to  survey  the  tumor  bed  and  assess  the  presence  of  residual  disease.  Last,  the  surface  of  the 
excised  specimen  is  imaged  with  radioluminescence  microscopy,  which  has  the  ability  to  visualize 
cancer  down  to  single  cells.  The  proposed  approach  will  therefore  provide  real-time  assessment  of 
the  disease  extent  at  different  stages  of  breast-conserving  surgery,  i.e.  before,  during,  and  after 
resection  of  the  primary  tumor.  This  multimodal  strategy  was  demonstrated  for  a  phantom  study  as 
wall  as  for  a  murine  breast  cancer  xenograft  model. 

Following  three  years  of  research  on  Cerenkov  luminescence,  we  can  draw  the  following  conclusions. 
Cerenkov  luminescence  is  a  practical  and  inexpensive  way  of  measuring  FDG  uptake  during  surgery. 
When  used  with  a  scintillator,  the  technique  allows  sensitive  detection  of  tracer  uptake,  down  to  a 
single  cell.  Thus  this  may  be  a  practical  strategy  for  precise  assessment  of  the  margins  of  disease. 
Sensitivity  is  limited  to  approximately  5  mm  depth,  which  means  that  Cerenkov  luminescence  imaging 
is  only  sensitive  to  superficial  malignancies,  which  is  an  advantage  for  assessing  the  status  of  the 
tissue  margins.  Last,  FDG  is  a  clinically  approved  radiotracer,  which  greatly  facilitates  translation  of 
this  proposed  imaging  technique. 

Although  Cerenkov  luminescence  imaging  is  a  promising  new  path  for  guiding  surgical  resection, 
there  are  two  main  issues  that  remain.  The  first  one  is  that  the  method  requires  the  administration  of 
a  radioactive  tracer  and  thus  it  exposes  the  patient  and  staff  to  low  level  ionizing  radiation.  Doctors 
and  patients  may  therefore  be  reluctant  to  undergo  such  procedures  (although  the  use  of  radiotracers 
is  standard  for  sentinel  lymph  node  biopsy).  The  other  issue  is  that  18F  decays  with  a  half-life  of  2h, 
thus  imaging  is  only  possible  within  a  2-3  hour  window,  which  could  be  problematic  for  some  of  the 
longer  surgeries. 
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Abstract 

Radiotracers  play  an  important  role  in  interrogating  molecular  processes  both  in  vitro  and  in  vivo.  However,  current  methods 
are  limited  to  measuring  average  radiotracer  uptake  in  large  cell  populations  and,  as  a  result,  lack  the  ability  to  quantify  cell- 
to-cell  variations.  Here  we  apply  a  new  technique,  termed  radioluminescence  microscopy ,  to  visualize  radiotracer  uptake  in 
single  living  cells,  in  a  standard  fluorescence  microscopy  environment.  In  this  technique,  live  cells  are  cultured  sparsely  on 
a  thin  scintillator  plate  and  incubated  with  a  radiotracer.  Light  produced  following  beta  decay  is  measured  using  a  highly 
sensitive  microscope.  Radioluminescence  microscopy  revealed  strong  heterogeneity  in  the  uptake  of  [18F]fluoro- 
deoxyglucose  (FDG)  in  single  cells,  which  was  found  consistent  with  fluorescence  imaging  of  a  glucose  analog.  We  also 
verified  that  dynamic  uptake  of  FDG  in  single  cells  followed  the  standard  two-tissue  compartmental  model.  Last,  we 
transfected  cells  with  a  fusion  PET/fluorescence  reporter  gene  and  found  that  uptake  of  FHBG  (a  PET  radiotracer  for 
transgene  expression)  coincided  with  expression  of  the  fluorescent  protein.  Together,  these  results  indicate  that 
radioluminescence  microscopy  can  visualize  radiotracer  uptake  with  single-cell  resolution,  which  may  find  a  use  in  the 
precise  characterization  of  radiotracers. 
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Introduction 

The  use  of  radiotracers  to  probe  biological  processes  has  several 
advantages  over  other  approaches:  radiotracers  can  be  synthesized 
with  chemical  composition  nearly  identical  as  a  given  compound 
of  interest;  their  concentration  measured  with  exquisite  sensitivity 
[1];  and  their  distribution  imaged  in  vivo  with  positron  emission 
tomography  (PET)  or  single  photon  emission  computer  tomogra¬ 
phy  (SPEGT)  [2, 3, 4, 5].  With  the  widespread  use  of  radionuclide 
imaging  in  research  and  in  hospitals,  we  need  to  better  understand 
how  properties  specific  to  individual  cells  (e.g.  gene  expression,  cell 
cycle,  cell  damage,  and  cell  morphology)  affect  the  uptake  and 
retention  of  radiotracers.  In  particular,  disease  and  therapy  can 
alter  cellular  mechanisms  in  a  heterogeneous  manner;  how  these 
alterations  affect  radiotracer  uptake  at  the  single-cell  level  is 
currently  unknown  and  of  critical  importance. 

The  averaging  effect  of  measuring  radiotracer  uptake  in  pooled 
cell  populations  can  mask  important  differences  between  cells 
belonging  to  the  same  population.  However,  current  approaches 
lack  the  ability  to  distinguish  radiotracer  uptake  in  individual  living 
cells.  For  instance,  film  autoradiographs  can  be  examined  with 
light  [6]  or  electron  microscopy  [7]  to  visualize  radioactive  decay 
within  individual  cells  but  the  method  is  limited  to  fixed  tissues  and 
low  energy  radionuclides  (e.g.  14C  and  3H).  Digital  autoradiog¬ 
raphy  techniques  (e.g.  storage  phosphor  [8],  electronic  detection 
[9,10],  thin  phosphor  layer  [11],  scintillator  [12],  and  gaseous 


chamber  [13])  offer  higher  detection  efficiency  and  dynamic  range 
but  poorer  spatial  resolution  (>30  flm),  insufficient  to  resolve 
individual  cells.  Likewise,  in  vivo  radiotracer  imaging  and 
scintillation  counting  can  only  measure  signals  from  large  cell 
populations. 

Here  a  new  method,  termed  radioluminescence  microscopy,  is 
proposed  to  measure  radiotracer  uptake  in  single  living  cells. 
Radioluminescence  is  the  physical  process  by  which  ionizing 
charged  particles  produce  light  in  certain  materials.  Due  to  the 
short  range  of  beta  particles  (electrons  or  positrons),  radiolumines¬ 
cence  occurs  near  the  location  of  the  radioactive  emitter.  The 
range  of  these  particles  is  further  reduced  in  dense,  high-atomic- 
number  materials  such  as  inorganic  scintillators.  Following  this 
observation,  we  hypothesized  that  the  radioactivity  of  single  cells 
could  be  measured  by  placing  these  cells  in  contact  with 
a  scintillator  plate  and  imaging  the  resulting  optical  signal  using 
a  sensitive  microscope  with  high  numerical  aperture  (NA)  and  high 
photon  sensitivity.  Furthermore,  we  envisioned  that  this  technique 
could  be  applied  concurrently  with  standard  fluorescence  micros¬ 
copy  because  scintillator  materials  are  optically  clear  in  the  visible 
range. 

The  proposed  radioluminescence  microscopy  set-up  consists  of 
a  100  Jim-thin  CdW04  scintillator  plate,  on  which  cells  have 
adhered,  immersed  in  a  glass-bottom  dish  filled  with  cell  culture 
medium  (Figure  1A&B).  The  dish  is  imaged  using  an  inverted 
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Figure  1.  Overview  of  the  radioluminescence  microscope.  (A)  Radioluminescence  is  produced  within  a  scintillator  plate  following  the 
emission  of  a  beta  particle  from  a  radiotracer  within  a  cell  (yellow  glow).  The  optical  photons  are  captured  by  a  high-numerical-aperture  objective 
coupled  to  a  deep-cooled  EM-CCD  camera.  Emission  and  excitation  filters  used  in  combination  with  a  light  source  allow  for  concurrent  fluorescence 
and  brightfield  microscopy.  (B)  Photograph  of  the  system  showing  a  glass-bottom  dish  containing  a  scintillator  plate  immersed  in  cell  culture 
medium  and  placed  into  the  inverted  microscope.  (C)  Three  GFP-expressing  HeLa  cells  located  near  the  corner  of  a  scintillator  plate  were  localized 
using  fluorescence  microscopy  (arrows).  The  edge  of  the  scintillator  plate  is  outlined  in  red.  (D)  After  incubation  with  FDG  (400  juCi,  1  h),  these  three 
cells  also  produced  focal  radioluminescence  signal  coincident  with  the  fluorescent  emission. 
doi:1 0.1 371 /journal. pone.0046285.g001 


microscope  fitted  with  a  high-NA  objective  and  an  electron- 
multiplying  charge-coupled  device  (EM-CCD). 

As  an  illustration  of  the  methods,  human  ovarian  cancer  cells 
(HeLa)  expressing  the  green  fluorescent  protein  (GFP)  were 
imaged  after  incubation  with  [18F]fluorodeoxyglucose  (FDG; 
400  pCi).  Three  isolated  cells  were  localized  near  the  corner  of 
a  scintillator  plate,  which  is  clearly  visible  on  the  brightfield 
micrograph  (Figure  1C,  dashed  red  line).  Both  fluorescence  and 
radioluminescence  images  displayed  focal  signal  at  the  locations  of 
the  three  cells  (arrows,  Figure  1C&D). 

Results 

Radioluminescence  Imaging  of  FDG  uptake  in  Single 
Cells 

FDG  is  preferentially  taken  up  and  retained  within  tissues  with 
high  glucose  metabolism  such  as  malignant  tumors  [14,15,16]. 
Measuring  FDG  uptake  in  a  heterogeneous  cell  population  is  of 
great  interest  as  it  may  help  better  understand  the  heterogeneous 
metabolic  alterations  displayed  by  tumors,  and  the  impact  that  the 
tumor  microenvironment  has  on  these  alterations  [17,18]. 
However,  there  does  not  exist  a  standard  method  for  measuring 
radiotracer  uptake  at  the  single  cell  level.  Therefore,  to  validate  the 
use  of  radioluminescence  microscopy  for  FDG  imaging,  we  used 
a  fluorescent  glucose  analog  as  a  surrogate  for  FDG  uptake  in 
single  cells. 

After  a  1  h  glucose  fasting  period,  we  incubated  human  breast 
cancer  cells  (MDA-MB-231)  for  1  h  at  37°C  with  FDG  (400  pCi) 
and  2-  [N-(7-nitrobenz-2-oxa- 1 ,3-diaxol-4-yl)amino]  -2-deoxyglu- 
cose  (2-NBDG;  100  pM)  [19,20].  After  washing  the  cells,  we 
acquired  brightfield,  radioluminescence  and  fluorescence  micro¬ 
graphs.  We  observed  good  co-localization  between  the  radiolumi¬ 
nescence  intensity  and  the  cell  outline  seen  on  brightfield  images 
(Figure  2A).  Furthermore,  the  radioluminescence  intensity  varied 
significantly  from  cell  to  cell,  indicating  heterogeneous  uptake  of 
FDG.  The  single-cell  radioluminescence  signal  was  correlated  with 


2-NBDG  fluorescence  (Figure  2B,  p<  10  5,  r=0.74).  An  exact 
correlation  between  FDG  and  2-NBDG  is  not  expected  due  to  (i) 
possibly  distinct  transport  me  chanisms  [2 1] ;  and  (ii)  the  inability 
of  2-NBDG  to  fluoresce  after  being  metabolized  [22].  A  line 
profile  through  the  fluorescence  and  radioluminescence  images 
confirms  co-localization  of  FDG  and  2-NBDG  signals 
(Figure  2C). 

Pharmacokinetic  Analysis  of  FDG  Metabolism  in  Single 
Cells 

The  transport  and  retention  of  FDG  in  a  cell  is  influenced  by 
multiple  factors,  such  as  the  expression  of  various  genes,  the 
density  of  glucose  transporters  on  the  cell  surface,  the  cell  size,  and 
the  levels  and  activities  of  hexokinase  and  phosphatase  enzymes 
[15].  Under  steady-state  conditions,  the  intracellular  and  extra¬ 
cellular  FDG  concentrations  are  in  equilibrium.  However,  rapid 
changes  in  the  extracellular  environment  induce  a  transient 
response  characteristic  of  the  cell’s  glucose  metabolism  parame¬ 
ters.  These  parameters  can  be  estimated  using  pharmacokinetic 
modeling  techniques.  The  ability  to  manipulate  a  cell’s  environ¬ 
ment  is  unique  to  an  in  vitro  setting  and  cannot  be  easily  replicated 
in  vivo.  Furthermore,  pharmacokinetic  modeling  from  PET  or 
gamma  counting  measurements  requires  assumptions  such  as 
uniform  radiotracer  concentration  and  homogeneous  rate  param¬ 
eters  for  each  compartment  [23].  These  assumptions  may  not  be 
satisfied  in  practice  because  each  cell  in  the  compartment  is 
characterized  by  unique  parameters.  Pharmacokinetic  modeling  at 
the  single-cell  level  may  provide  more  optimal  characterization  of 
cellular  parameters. 

To  investigate  the  utility  of  radioluminescence  microscopy  for 
single-cell  pharmacokinetic  studies,  we  monitored  the  uptake  of 
FDG  in  breast  cancer  cells  (MDA-MB-231)  over  8  h.  After 
depriving  cells  of  glucose  for  1  h,  we  added  FDG  (5  pCi)  to  their 
medium  and  acquired  serial  brightfield  and  radioluminescence 
images  every  6  min  for  8  h  (Figure  3A  &  Video  SI).  Although 
FDG  uptake  varied  significantly  from  cell  to  cell,  all  cells  displayed 
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Figure  2.  Radioluminescence  imaging  of  FDG  uptake  in  single  cells.  Human  breast  cancer  cells  (MDA-MB-231)  were  deprived  of  glucose  for 
1  h,  incubated  for  1  h  with  FDG  (400  pCi)  and  2-NBDG  (100  juM),  and  then  washed.  (A)  Brightfield  (scale  bar,  100  pm.),  radioluminescence  (FDG),  and 
fluorescence  (2-NBDG)  micrographs  (Objective:  40X/1.3  NA).  Overlay,  showing  co-localized  radioluminescence  (green)  and  fluorescence  (red).  (B) 
Scatter  plot  comparing  FDG  and  2-NBDG  uptake,  computed  over  140  cells  (light  red  dots)  and  26  control  ROIs  (blue  dots).  The  green  line  was 


obtained  by  linear  regression  (correlation,  r  =  0.74).  Arbitrary  units  (A.U.). 
along  a  line  profile  [red  dashed  line  in  (A)]. 
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the  same  linear  increase  in  radioactivity,  followed  by  a  plateau  and 
a  slow  decrease  after  3  h  (Figure  3D). 

We  performed  two  other  sets  of  experiment  to  highlight  efflux  of 
FDG  from  a  cell.  Toward  this  goal,  we  next  subjected  breast 
cancer  cells  (MDA-MB-231)  to  conditions  known  to  minimize 
FDG  influx,  i.e.  competition  from  glucose  (Figure  3B  &  Video 
S2)  and  withdrawal  of  FDG  (Figure  3C  &  Video  S3).  The 
addition  of  glucose  to  the  medium  (25  mM)  at  2  h  lead  to  a  strong 
decline  in  cell  radioactivity  (Figure  3E)  as  FDG  and  glucose 
competed  for  the  same  glucose  transporters.  Withdrawing  FDG 
from  the  media  of  cells  that  had  previously  been  incubated  with 
FDG  (400  pCi,  1  h)  also  resulted  in  a  similarly  fast  decrease  in  cell 
radioactivity  (Figure  3F). 

The  uptake  and  metabolism  of  FDG  can  be  mathematically 
modeled  using  a  two-tissue  compartmental  model  (Figure  4A), 
whose  rate  constants  K\,  Jc2,  k 3  and  £4  represent  the  influx,  efflux, 
phosphorylation,  and  dephosphorylation  of  FDG,  respectively. 
Influx  of  FDG  in  cells  (as  shown  in  Figure  3A)  was  quantified  by 
Patlak  analysis.  Single-cell  time-activity  curves  measured  by 
radioluminescence  microscopy  were  found  consistent  with  Patlak’s 
model,  at  least  in  the  early  time  points:  After  a  short  transient 
period,  equilibrium  was  established  and  the  intracellular  concen¬ 
tration  of  FDG  increased  linearly  with  time  due  to  the  irreversible 
trapping  of  FDG  into  the  cell  (e.g.  Figure  4B).  The  slope  of  the 
linear  rise  is  the  product  of  two  terms,  namely  K\ ,  the  influx  rate, 
and  £3/(^2  +  ^3)5  the  fraction  of  the  intracellular  FDG  irreversibly 
metabolized. 


PLOS  ONE  |  www.plosone.org 


(C)  Radioluminescence  (FDG)  and  fluorescence  (2-NBDG)  intensity  shown 


We  found  large  variations  in  the  Patlak  coefficients  across  the 
cells  that  were  imaged,  indicating  that  seemingly  identical  cells 
process  glucose  heterogeneously.  Furthermore,  solving  for  the 
pharmacokinetic  coefficients  K\ ,  k2  and  k 3  showed  that  K\  (influx) 
and  k2  (efflux)  were  correlated  (/?  <  10-5,  r  =  0.89,  Figure  4C)  but 
K\  and  k 3  (phosphorylation)  were  not  (p  =  0.6,  r=  —  0.08, 
Figure  4D).  Also,  the  majority  of  cells  stopped  accumulating 
FDG  at  approximately  3  h  and  a  slow  decrease  in  cell  FDG 
concentration  was  observed  (Figure  3D).  The  non-negligible  rate 
of  FDG  dephosphorylation  (Aq)  is  likely  the  main  factor 
contributing  to  that  effect.  However,  dephosphorylation  alone 
should  result  in  the  FDG  concentration  reaching  a  steady  plateau 
due  to  equilibration  of  phosphorylation  and  dephosphorylation. 
The  slow  decrease  that  was  observed  instead  may  have  been 
caused  by  increased  competition  from  unlabeled  2DG  (a 
byproduct  of  FDG  synthesis)  as  FDG  concentration  diminished 
due  to  radioactive  decay. 

We  also  derived  a  mathematical  model  to  represent  FDG  efflux 
from  a  cell  after  withdrawal  of  FDG  (as  shown  in  Figure  3C), 
composed  of  the  sum  of  a  slow  and  a  fast  exponential  decay.  The 
model  was  found  to  be  in  agreement  with  radioluminescence 
measurements  of  single  cells  (e.g.  Figure  4E),  confirming  that  two 
processes  are  occurring  concurrently  at  different  rates.  The  first 
process  describes  the  rapid  diffusion  of  unbound  FDG  out  of  the 
cell  (rate  fa),  whereas  the  second  process  involves  the  slow 
dephosphorylation  of  FDG-6-phosphate  (rate  fa).  While  the  efflux 
rate  was  heterogeneous  over  the  cell  population  studied,  we  found 
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Figure  3.  Dynamic  radioluminescence  imaging  of  FDG  in  single  cells.  Micrographs  (brightfield  and  radioluminescence)  were  acquired  every 
6  min  for  8  h  for  three  experiments.  (A)  MDA-MB-231  cells  are  imaged  while  being  incubated  with  FDG  (5  juCi).  (B)  Glucose  (25  mM)  is  added  2  h 
after  the  beginning  of  the  incubation  with  FDG  (5  jiCi).  (C)  FDG  is  withdrawn  at  the  start  of  imaging  after  incubation  (1  h,  400  jLiCi).  Scale  bar:  100  pm. 
(E-F)  Time-activity  curves  plotted  for  individual  cells  (light  red  lines)  and  1 0  control  ROIs  manually  selected  in  the  background  (light  blue  lines),  for  all 
three  experiments.  The  thick  red  and  blue  lines  represent  the  average  for  cells  and  control  ROIs,  respectively. 
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no  significant  correlation  between  the  fast  and  slow  components  of 
the  decay  (p  =  0.9,  r=  —0.02,  Figure  4F). 

Single-cell  Characterization  of  Transgene  Expression  with 
a  PET  Probe 

To  further  validate  radioluminescence  microscopy,  we  in¬ 
vestigated  the  uptake  of  9-(4-[18F]Fluoro-3-hydroxymethylbutyl)- 
guanine  (FHBG)  in  cancer  cells  that  were  heterogeneously 
transfected  to  express  the  mutant  herpes  simplex  virus  type  1 
truncated  thymidine  kinase  (HSVl-ttk).  HSVl-ttk  can  selectively 
metabolize  and  trap  radiolabeled  substrates  such  as  FHBG  [24]. 
Because  FHBG  has  low  affinity  for  mammalian  thymidine  kinases 
(TK)  and  high  affinity  for  viral  HSV1-TK,  it  can  be  used  to  image 
cell  trafficking  in  living  subjects  with  PET  [25].  To  assess  the 


expression  of  the  HSVl-ttk  transgene  with  fluorescence  micros¬ 
copy,  we  built  a  fusion  reporter  that  also  encodes  the  monomeric 
red  fluorescent  protein  1  (mrfpl). 

We  transfected  human  cervical  cancer  cells  (HeLa)  with  the 
fusion  reporter  vector  encoding  HSVl-ttk  and  mrfpl.  Radiolu¬ 
minescence  microscopy  of  FHBG  (incubation  2h  with  300  pCi) 
demonstrated  focal  radiotracer  uptake,  with  individual  cells  clearly 
resolvable  under  100X  magnification  (Figure  5A).  Using  fluores¬ 
cence  microscopy,  we  estimated  that  88%  of  the  cells  (217/245) 
had  been  successfully  transfected  with  the  fusion  transgene.  All  of 
those  cells  were  also  clearly  distinguishable  on  the  radiolumines¬ 
cence  fluorescence  image  (Figure  5B).  We  also  found  that  9%  of 
the  cells  (21/245)  had  not  been  transfected  and  did  not  produce 
any  fluorescence.  These  cells  did  not  present  a  radioluminescence 
signal,  which  suggest  that  they  did  not  retain  FHBG.  The 
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Figure  4.  Pharmacokinetics  analysis  in  single  cells.  (A)  Two-tissue  compartmental  model  describing  FDG  pharmacokinetics,  including  influx 
{K i),  efflux  (k2),  phosphorylation  to  FDG-6-phosphate  (k3),  and  dephosphorylation  (/c4).  (B)  Patlak  analysis  modeling  FDG  influx  kinetics  for  a  single  cell 
(highlighted  by  a  red  circle  in  Figure  3A).  (C,D)  Rate  of  efflux  (/c2)  and  phosphorylation  (/c3)  plotted  as  a  function  of  rate  of  influx  (/Cn)  for  all  the  cells  in 
the  microscope's  field  of  view.  (E)  Compartmental  analysis  modeling  FDG  efflux  kinetics  from  a  single  cell  (highlighted  by  a  red  circle  in  Figure  3C) 
after  withdrawal  of  FDG,  presenting  a  fast  and  a  slow  component.  (F)  The  model  for  FDG  efflux  is  the  sum  of  a  fast  and  a  slow  component  (rates  2-| 
and  /2,  respectively),  which  are  plotted  for  all  the  cells  in  the  field  of  view. 
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remaining  5%  of  the  cells  (7/245)  were  excluded  from  the  analysis 
due  to  ambiguous  radioluminescence  intensity,  mostly  due  to  the 
proximity  of  one  or  more  strongly  positive  cells  (e.g.  Figure  5B, 
green  arrow).  Generally,  in  our  system,  radioluminescence  signals 
for  FHGB-positive  and  negative  cells  were  more  distinctly 
separated  than  fluorescence  signals  for  RFP-positive  and  negative 
cells  (Figure  5B,  white  arrows).  A  line  profile  passing  through  four 
cells  showed  good  co-localization  of  RFP  and  FHBG  (Figure  5D). 

While  uptake  of  FHBG  was  coincident  with  RFP  fluorescence, 
fluorescence  intensity  was  not  strongly  predictive  of  radiolumines¬ 
cence  intensity  (Figure  5C;  correlation,  r  =  0.34),  indicating  that 
although  the  HSVl-tk  reporter  gene  expression  is  required  for 
FHBG  uptake,  the  level  of  transgene  expression  is  not  solely 
responsible  for  the  extent  of  FHBG  uptake.  In  a  separate 
experiment,  the  FHBG  substrate  displayed  no  affinity  for 
mammalian  TK  enzyme:  wild-type  HeLa  cells  incubated  with 
FHBG  (300  pCi,  2h)  showed  no  measureable  radioluminescence 
signal  (Figure  5E). 

Performance  Characterization 

To  investigate  the  spatial  resolution  of  the  imaging  set-up,  dry 
FDG  aggregates  were  imaged  with  the  radioluminescence 
microscope.  Brightfield  and  radioluminescence  images  displayed 
good  correlation  (r=  —  0.79,  /?<10-5;  Figure  6A-G).  From 
these  measurements,  we  estimated  the  microscope  spatial  resolu¬ 
tion  to  be  5  pm  (full-width  half-maximum). 

The  sensitivity  of  the  microscope  was  measured  by  imaging  the 
decay  of  a  uniform  distribution  of  FDG  (2.6  pCi  initially)  over  24 
hours.  The  average  pixel  intensity  (expressed  as  a  percentage  of 
the  pixel  intensity  in  the  first  frame)  decreased  exponentially  with 
time  with  a  half-life  consistent  with  the  decay  of  18F  (Figure  6H). 


The  relationship  between  the  average  pixel  intensity  and  the 
average  activity  per  area  remained  linear  all  the  way  down  to 
approximately  0.1  fCi/pm2.  The  signal-to-noise  ratio  (SNR)  also 
decreased  with  decreasing  activity  (Figure  61).  To  visualize  image 
features,  a  SNR  of  at  least  five  is  required  (Rose  criterion  [26]), 
which  corresponds  to  a  minimum  activity  area  density  of  4  fCi/ 
pm2.  This  is  equivalent  to  1.4  molecules  of  FDG  per  pm2. 

Discussion 

For  the  first  time,  radioluminescence  microscopy  enables  the 
quantification  of  radiotracer  uptake  and  pharmacokinetics  at  the 
single-cell  level.  Because  a  similar  method  does  not  currently  exist, 
we  validated  the  new  approach  by  comparing  single-cell  radio- 
tracer  measurements  against  surrogate  quantities.  Hence,  we 
found  FDG  uptake  to  be  consistent  with  fluorescence  imaging  of 
a  glucose  analog  (Figure  2),  with  pharmacokinetics  characteristic 
of  a  two-tissue  compartmental  model  (Figure  3  and  Figure  4). 
Furthermore,  in  a  cell  population  heterogeneously  transfected  with 
a  fusion  PET/fluorescence  reporter,  we  verified  that  FHBG 
uptake  was  concordant  with  fluorescence  imaging  of  mrfp  1 

(Fig-  5). 

Although  radioluminescence  microscopy  is  not  yet  capable  of 
visualizing  intracellular  radiotracer  distributions,  it  can  measure 
the  radioactivity  of  single  cells  provided  that  those  cells  are 
spatially  separated  on  the  scintillator  plate.  Accurate  measure¬ 
ments  can  be  achieved  with  cell-to-cell  separation  of  10  pm  or 
more.  Radioluminescence  micrographs  can  be  acquired  in  5  min 
or  less  using  mostly  off-the-shelf  instrumentation.  Commonly  used 
beta  emitters  such  as  18F,  131I,  and  64Cu  can  be  used  to  produce 
such  images. 
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Figure  5.  Radioluminescence  imaging  of  gene  expression  in  single  cells.  Human  cervical  cancer  cells  (HeLa)  transfected  with  a  fusion  PET/ 
fluorescence  reporter  gene  were  incubated  with  FHBG  (300  pCi,  2  h).  (A)  Brightfield  (scale  bar,  50  pm),  radioluminescence  (FHBG),  and  fluorescence 
(RFP)  micrographs  (objective,  100X/1.35  NA).  Overlay  shows  FHBG  radioluminescence  (green),  RFP  fluorescence  (red),  and  cell  outline  segmented 
from  brightfield.  Cells  negative  for  RFP  are  also  negative  for  FHBG  (red  arrows).  (B)  Same  as  (A),  but  with  a  40X/1.3  NA  objective  (scale  bar,  100  pm). 
White  arrows  indicate  cells  with  weak  fluorescence  intensity  but  substantial  radioluminescence  intensity.  The  green  arrow  points  to  a  cell  with  no  RFP 
expression  but  ambiguous  radioluminescence  intensity.  (C)  Scatter  plot  of  FHBG  vs.  RFP  uptake,  computed  for  245  cells  (light  red  dots)  and  100 
control  ROIs  (blue  dots).  Arbitrary  units.  (D)  Radioluminescence  and  fluorescence  shown  along  a  line  profile  [red  dashed  line  in  (A)].  (E)  Same 
experiment  as  (A,B),  but  using  control  wild-type  HeLa  cells  (scale  bar,  100  pm). 
doi:1 0.1 371 /journal. pone.0046285.g005 


While  radioluminescence  microscopy  is  mainly  intended  to 
image  tissue  culture  cells,  the  method  may  be  applicable  to 
imaging  solid  tissue  sections.  However,  in  the  current  configura¬ 
tion,  it  may  not  provide  single-cell  resolution  for  dense  tissue 
section.  One  solution  to  this  problem  is  to  dissociate  the  tissue 
prior  to  imaging  to  ensure  sufficient  separation  between  cells  [27]. 
We  are  also  currently  investigating  several  approaches  to  further 
improve  the  spatial  resolution  of  the  system. 

We  expect  that  radioluminescence  microscopy  will  become 
a  useful  technique  for  the  precise  characterization  of  radiotracer 
uptake  and  pharmacokinetics  at  the  single-cell  level.  New 
developments  in  scintillator  research  will  undoubtedly  improve 
the  performance  of  the  technique.  Thinner  scintillator  plates 
with  higher  density  and  light  yield  will  provide  better  spatial 
resolution  and  signal-to-noise  ratio.  Progress  in  image  processing 
and  calibration  techniques  will  also  allow  for  more  quantitative 
measurements  of  radiotracer  concentration  in  single  cells. 

Materials  and  Methods 

Microscopy  Set-up 

Adherent  cancer  cells  were  seeded  sparsely  on  5  mm  x  5  mm  x 
0.1  mm  plates  made  of  CdW04,  a  non-hygroscopic  inorganic 
scintillator,  with  both  sides  polished  (Figure  1A).  CdW04  has 
relatively  high  light  yield  (12,000-15,000  photon/MeV),  high 
effective  atomic  number  (Zeff=64),  high  density  (7.9  g/cm3),  and 
no  significant  afterglow.  The  scintillator  plates,  loaded  with  cells, 


were  placed  in  microscopy  dishes  (#0  cover  glass,  0.085- 
0.115  mm,  In  Vitro  Scientific)  filled  with  fresh  media 
(Figure  1A&B).  The  use  of  thin  scintillator  plates  and  thin- 
bottom  imaging  dishes  is  required  to  accommodate  the  short 
working  distance  of  the  microscope  objective  (200  pm). 

The  imaging  dishes  were  placed  in  a  bioluminescence  micro¬ 
scope  (LV200,  Olympus)  outfitted  with  either  a  40X/1.3  NA  oil 
objective  (UPLFLN40XO,  Olympus)  or  a  100X/1.35  NA  oil 
objective  (UPLAPOOOXOI3,  Olympus),  and  a  deep-cooled 
electron-multiplying  charge-coupled  device  (EM-GGD;  ImageEM 
C9100-14,  Hamamatsu)  (Figure  1  A&B).  The  C9100-14  GCD  is 
a  back-thinned  frame  transfer  device,  with  a  1024x1024  array  of 
13  pmxl3  pm  pixels.  The  LV200  is  also  equipped  with 
temperature,  humidity,  and  C02  regulation  for  extended  live  cell 
imaging. 

Brightfield  images  were  acquired  with  no  EM  gain,  a  neutral- 
density  filter  on  the  excitation,  and  the  emission  shutter  open.  For 
the  40X  magnification,  radioluminescence  images  were  taken  with 
an  exposure  time  of  5  min,  an  EM  gain  of  251/1200,  2x2  pixel 
binning,  the  excitation  shutter  closed,  and  the  emission  shutter 
open.  For  the  100X  magnification,  the  exposure  time  was  20  min 
and  the  EM  gain  505/1200.  We  used  the  brightfield  mode  to  set 
the  microscope  into  focus.  Optimal  radioluminescence  focus  was 
achieved  when  the  cells  displayed  sharp  positive  contrast  in  the 
corresponding  brightfield  image.  For  fluorescence  microscopy,  we 
used  a  460  nm/535  nm  filter  set  for  2-NBDG  imaging  (Chroma 


PLOS  ONE  |  www.plosone.org 


6 


October  2012  |  Volume  7  |  Issue  10  |  e46285 


Radioluminescence  Microscopy  of  Single  Cells 


Figure  6.  Performance  characterization.  FDG  aggregates  were  obtained  by  evaporating  an  aqueous  solution  of  FDG  between  a  scintillator  and 
a  glass-bottom  imaging  dish.  (A)  Fused  radioluminescence  and  brightfield  images;  (B)  Brightfield  and  (C)  radioluminescence  images,  magnified;  (D) 
Brightfield  and  (E)  radioluminescence  images,  further  magnified,  focusing  on  one  particular  FDG  aggregate;  (F,G)  2-D  Gaussian  fit  of  (D)  and  (E), 
respectively.  (H)  Radioluminescence  microscope  sensitivity,  obtained  by  imaging  the  decay  of  a  drop  of  FDG  (2.6  jLiCi)  over  time.  Solid  line :  mean  pixel 
intensity;  Dashed  line :  ideal  exponential  decay  for  18F.  (I)  Per-pixel  signal-to-noise  ratio,  defined  as  the  ratio  of  the  average  pixel  intensity  to  the  noise 
standard  deviation.  The  sensitivity  of  the  system  is  defined  here  as  the  amount  of  activity  required  per  area  to  achieve  a  SNR  of  5  (Rose  criterion  ). 
doi:1 0.1 371/journal,  pone.0046285.g006 


Technology  Corp.,  filter  ref.  D460/50x  and  D535/40  m)  and 
a  540  nm/600  nm  filter  set  for  RFP  imaging  (Chroma  Technol¬ 
ogy  Corp.,  filter  ref.  HQ540/40x  and  HQ600/50  m). 

Samples  Preparation 

MDA-MB-231  human  breast  cancer  cells  were  purchased  from 
the  American  Type  Culture  Collection  (ATCC,  Manassas,  VA) 
and  cultured  in  Leibovitz’s  LI 5  medium  supplemented  with  10% 
fetal  bovine  serum.  One  side  of  the  scintillator  plate  was  coated 
with  fibronectin  (10  pg/ml)  to  allow  the  cells  to  attach.  After  the 
plate  had  dried,  the  cells  were  seeded  by  placing  a  50  pi  drop 
containing  104  cells  on  the  fibronectin-coated  plate.  Imaging  was 
performed  the  following  day. 

PCR  amplification  and  standard  cloning  techniques  were  used 
to  insert  the  mrfp  and  ttk  genes  from  plasmid  pCDNA3.1-CMV- 
hrl-mrfp-ttk.  A  lentiviral  EFl-gfp  vector  was  purchased  from 
System  Biosciences  (SBI,  Mountain  view,  CA).  The  gfp  fragment 
was  removed  from  the  vector  and  replaced  with  mrfp-ttk.  For 
PCR  amplification,  different  5'  and  3'  end  primers  were  used  to 
generate  the  fusion  vector  (EF1 -mrfp-ttk). 

HeFa  human  cervical  cancer  and  293T  human  embryonic 
kidney  cells  were  purchased  from  ATCC  and  cultured  in  high- 
glucose  Dulbecco’s  modified  eagle  medium  supplemented  with 
10%  fetal  bovine  serum.  293T  cells  were  used  to  produce  the 
lentivirus  following  standard  procedures.  HeFa  cells  were 
transfected  with  concentrated  lentivirus  for  48  h,  then  trypsinized 
and  seeded  onto  a  scintillator  plate  coated  with  fibronectin  (10  pg/ 
ml)  one  day  before  imaging. 


Imaging  Protocol 

For  static  imaging  of  glucose  metabolism  using  combined 
fluorescence  and  radioluminescence  microscopy,  MDA-MB-231 
cells  were  deprived  of  glucose  in  Feibovitz’s  F-15  medium  1  h 
before  incubation  with  FDG  (400  pCi)  and  2-NBDG  (Invitro- 
gen,  100  pM).  FDG  was  produced  at  the  Stanford  radiochem¬ 
istry  facility  using  an  on-site  cyclotron.  Experiments  were 
conducted  shortly  after  synthesis  of  FDG  to  achieve  high 
specific  activity. 

For  dynamic  imaging  of  glucose  metabolism,  three  experi¬ 
ments  were  conducted:  In  the  first  experiment  (“FDG”),  MDA- 
MB-231  cells  were  deprived  of  glucose  for  1  h,  after  which 
5  pCi  of  FDG  was  added  to  their  medium.  Imaging  started 
a  few  minutes  later.  In  the  second  experiment  (“FDG+glucose”), 
the  same  procedures  were  followed.  Additionally,  25  mM  of 
glucose  was  added  to  the  medium  at  2  h.  In  the  third 
experiment  (“FDG  withdrawal”),  the  cells  were  preliminarily 
incubated  in  FDG  (400  pCi,  1  h),  and  imaging  started 
approximately  15  minutes  after  cell  washing. 

For  imaging  of  transgene  expression  with  FHBG  and  RFP, 
transfected  HeFa  cells  were  incubated  for  2  h  with  300  pCi  of 
18F-FHBG.  The  FHBG  substrate  was  produced  at  the  Stanford 
radiochemistry  facility  using  an  on-site  cyclotron.  To  gather 
sufficient  cell  numbers,  five  fields  were  imaged  at  40X  using 
radioluminescence,  fluorescence  and  brightfield  microscopy.  One 
field  was  also  imaged  with  a  100X  objective. 

Image  Corrections  and  Analysis 

Image  correction  and  analysis  were  performed  using  MATFAB 
R2010a  (Mathworks,  Natick,  MA).  Radioluminescence  micro¬ 
graphs  were  corrected  by  subtraction  of  a  dark  image,  taken  with 
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the  same  exposure  time  but  with  a  non-radioactive  sample  in  the 
microscope.  These  images  were  further  corrected  for  field  flatness 
using  a  flat-field  calibration  map  acquired  using  a  uniform 
distribution  of  FDG.  Gaussian  filtering  was  applied  where 
appropriate  to  reduce  noise.  During  long  exposures,  high-energy 
photons  (gamma  rays  and  annihilation  photons)  interacted  with 
the  CCD  and  produced  hot  spots  in  the  image.  These  hot  spots 
were  removed  by  applying  a  custom  algorithm  that  can  detect 
sharp  features  well  above  neighboring  pixels.  All  radiolumines¬ 
cence  images  were  corrected  for  radioactive  decay.  The  timestamp 
of  the  first  acquired  image  was  used  as  the  reference  time  point. 

Fluorescence  micrographs  were  corrected  for  background 
effects  (filter  bleed-through  and  camera  dark  noise)  by  subtracting 
a  dark  image  taken  with  a  non-fluorescent  sample.  Field  flatness 
was  corrected  using  a  flat-field  calibration  map. 

To  measure  radiotracer  uptake  in  single  cells,  circular  regions  of 
interest  (ROIs;  diameter,  24  pm)  were  manually  placed  on  the 
cells  using  the  brightfield  micrograph.  Similar  ROIs  were  placed 
in  the  background  as  controls.  Cell  radiotracer  uptake  was  defined 
as  the  mean  pixel  intensity  within  the  ROI  of  the  corrected 
radioluminescence  image.  The  same  ROI  analysis  procedure  was 
also  applied  to  fluorescence  micrographs. 

The  range  of  cell  motion  occurring  during  the  exposure  of 
a  single  frame  is  typically  too  small  to  result  in  any  significant 
blurring.  However,  for  extended  timelapse  imaging  studies  (1  h  or 
longer),  cell  motion  can  no  longer  be  neglected  and  must  be 
accounted  for  when  analyzing  the  images.  We  therefore  manually 
placed  circular  ROIs  on  the  each  cell  every  10  frames  (i.e.  every 
hour).  In  between  these  key  frames,  we  assumed  that  the  cells 
moved  in  a  straight  line. 

Radiotracer  Kinetic  Modeling 

Influx  of  FDG  into  glucose-deprived  cells  was  described  using 
the  following  two-tissue  compartmental  model: 


Ca 


Kiki 

(k2+k2)2 


(\-e-(h+k3  )*)  + 


Kik3  t 
k2  +  k  3 


where  Ca  is  the  extracellular  FDG  concentration  (assumed  to  be 
fixed);  C{t)  is  the  time-dependent  intracellular  FDG  concentration 
(including  free  FDG  and  bound  FDG-6-P);  and  K\,  k2  and  k2  are 
the  rate  constants  representing  influx,  efflux,  and  irreversible 
phosphorylation  of  FDG,  respectively  [23].  For  t>  >\/(k2-\-k$\ 
the  exponential  term  is  becomes  negligible.  The  intracellular  and 
extracellular  compartments  are  then  in  equilibrium,  with  the 
intracellular  concentration  of  FDG  rising  linearly  with  time  due  to 
irreversible  trapping.  The  slope  and  intercept  of  this  linear  rise  are 
the  Patlak  coefficients  [28].  We  used  non-linear  weighted  least- 
squares  curve  fitting  to  estimate  the  parameters  of  the  model.  The 
fitting  weights  were  adjusted  to  decrease  the  contribution  of  later 
time  points,  which  have  higher  noise  due  to  radioactive  decay. 

Efflux  of  FDG  from  cells  was  modeled  using  a  two-tissue 
compartmental  model: 


C(t)  =  a\e  Al t  +  a2e  Xlt 

where  a\  and  a2  are  positive  coefficients  that  depend  on  the  initial 
conditions,  and  k\  and  k2  are  the  eigenvalues  of  the  differential 
system  of  equations  describing  transport  of  FDG  between 
compartments.  The  rate  constant  Aq,  which  models  the  possible 
dephosphorylation  of  FDG-6-phosphate  (FDG-6-P),  was  included 
in  this  model  but  assumed  to  be  much  smaller  than  k2. 
Furthermore,  due  to  the  large  extracellular  volume  (0.2  ml),  the 


concentration  of  FDG  in  the  cell  culture  medium  was  assumed  to 
remain  negligible  after  withdrawal  of  FDG.  Under  these 
assumptions,  the  eigenvalues  can  be  approximated  as 


k\  —  k2-\-k2-\- 


k2k/[ 
k2  +  k2 


and 


_  k2k4 
2  k2  +  k 3- 


These  rate  parameters  were  estimated  by  fitting  the  efflux  model 
to  the  measured  time-activity  curves.  For  cells  for  which  the 
solution  of  the  fit  yielded  k\  ~k2  or  k2  <  1  min-1,  the  efflux  curve 
was  fitted  with  a  single  exponential  function.  In  the  special  case  of 
irreversible  trapping  (Aq  =  0),  the  model  is  described  by  a  single 
exponential  decay  with  rate  k\=k2+k2. 

Statistical  Analysis 

Correlation  between  fluorescence  and  radioluminescence  ROI 
measurements  was  computed  using  the  Pearson  product-moment 
correlation  coefficient.  A  p-value  of  less  than  0.01  was  considered 
statistically  significant. 

Spatial  Resolution  Characterization 

To  evaluate  the  performance  of  the  radioluminescence  micro¬ 
scope,  we  placed  a  drop  of  FDG  (10  pCi)  between  the  imaging 
dish  and  a  scintillator  plate.  Upon  evaporation  of  the  aqueous 
solvent,  FDG  precipitated  into  small  solid  aggregates  that  could  be 
seen  both  on  brightfield  and  radioluminescence  images.  We 
measured  the  size  of  these  aggregates  by  fitting  them  with  2-D 
Gaussian  functions. 

Sensitivity  Characterization 

The  overall  sensitivity  of  the  system  was  evaluated  by  imaging 
the  decay  of  a  mixture  of  glycerol  and  FDG  (2.6  pCi  initially), 
placed  between  an  imaging  dish  and  a  scintillator  plate. 
Radioluminescence  images  were  acquired  every  31  min,  using 
an  EM  gain  of  251/1200  and  an  exposure  time  of  30  min.  Within 
a  large  region  of  interest  (370,000  pixels),  pixel  intensities  were 
normalized  to  correct  for  field  flatness  using  the  first  frame  as 
a  reference.  The  standard  deviation  of  the  noise,  found  to  be 
Gaussian-distributed,  was  computed  in  each  frame.  The  per-pixel 
signal  to  noise  ratio  was  then  defined  as  the  ratio  of  the  average 
pixel  intensity  to  the  standard  deviation  of  the  noise.  The 
sensitivity  of  the  system  was  defined  as  the  activity  (area  density) 
required  to  achieve  a  signal-to-noise  ratio  of  5  (Rose  criterion 
[26]). 

Supporting  Information 

Video  SI  Timelapse  imaging  of  FDG  influx  kinetics  in 
MDA-MB-231  cells  using  radioluminescence  {left)  and 
brightfield  {right)  microscopy.  The  cells  were  deprived  of 
glucose  one  hour  prior  to  imaging.  Serial  image  acquisition  was 
started  after  adding  FDG  (5  pCi)  to  the  cell  culture  medium. 
(MP4) 

Video  S2  Timelapse  imaging  of  competition  between 
FDG  and  glucose  uptake  in  MDA-MB-231  cells  using 


PLOS  ONE  |  www.plosone.org 


8 


October  2012  |  Volume  7  |  Issue  10  |  e46285 


Radioluminescence  Microscopy  of  Single  Cells 


radioluminescence  ( left )  and  brightfield  [right)  micros¬ 
copy.  The  cells  were  deprived  of  glucose  one  hour  prior  to 
imaging.  Serial  image  acquisition  was  started  after  adding  FDG 
(5  jlCi)  to  the  cell  culture  medium.  Additionally,  25  mM  of 
glucose  was  added  to  the  medium  at  2  h. 

(MP4) 

Video  S3  Timelapse  imaging  of  FDG  efflux  kinetics  in 
MDA-MB-23 1  cells  using  radioluminescence  [left)  and 
brightfield  [right)  microscopy.  The  cells  were  preliminarily 
incubated  in  FDG  (400  JdCi,  1  h),  and  imaging  started  approxi¬ 
mately  15  minutes  after  cell  washing. 

(MP4) 
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Abstract.  Monte  Carlo  simulation  is  considered  the  most  reliable  method  for  modeling  photon  migration  in 
heterogeneous  media.  However,  its  widespread  use  is  hindered  by  the  high  computational  cost.  The  purpose  of 
this  work  is  to  report  on  our  implementation  of  a  simple  MapReduce  method  for  performing  fault-tolerant  Monte 
Carlo  computations  in  a  massively-parallel  cloud  computing  environment.  We  ported  the  MC321  Monte  Carlo 
package  to  Hadoop,  an  open-source  MapReduce  framework.  In  this  implementation,  Map  tasks  compute  photon 
histories  in  parallel  while  a  Reduce  task  scores  photon  absorption.  The  distributed  implementation  was  evaluated 
on  a  commercial  compute  cloud.  The  simulation  time  was  found  to  be  linearly  dependent  on  the  number  of 
photons  and  inversely  proportional  to  the  number  of  nodes.  For  a  cluster  size  of  240  nodes,  the  simulation  of  1 00 
billion  photon  histories  took  22  min,  a  1258  x  speed-up  compared  to  the  single-threaded  Monte  Carlo  program. 
The  overall  computational  throughput  was  85,1  78  photon  histories  per  node  per  second,  with  a  latency  of  1 00  s. 
The  distributed  simulation  produced  the  same  output  as  the  original  implementation  and  was  resilient  to  hardware 
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1  Introduction 

Researchers  have  long  relied  on  single-threaded  programming 
for  solving  computational  problems.  However,  in  many  appli¬ 
cations,  the  growth  of  scientific  data  has  outpaced  the  perfor¬ 
mance  of  single-core  processors.  Furthermore,  multicore  pro¬ 
cessors  and  many-core  graphics  processing  units  (GPUs)  are 
now  the  industry  standard  for  high-performance  computing.  As 
a  result,  scientific  computing  is  inexorably  shifting  to  parallel 
architectures. 

In  biomedical  applications,  the  need  for  high-performance 
computing  is  growing  as  technology  evolves  toward  more  ac¬ 
curate  imaging  and  treatment  delivery  methods.  Monte  Carlo 
simulation  is  the  gold  standard  for  modeling  complex  physi¬ 
cal  systems,  such  as  photon  migration  in  biological  tissue.1-3 
Yet  the  use  of  Monte  Carlo  methods  is  still  hampered  by  the 
high  computational  cost.  For  instance,  it  is  still  not  practical 
to  conduct  Monte  Carlo  simulations  for  planning  of  photody¬ 
namic  therapy,4  especially  when  using  inhomogeneous  tissue 
models.5,6 

Several  distributed  approaches  have  been  proposed  and  im¬ 
plemented  to  accelerate  Monte  Carlo  simulation  of  photon  trans¬ 
port.  Monte  Carlo  simulation  belongs  to  a  class  of  problems 
referred  to  as  “embarrassingly  parallel”,  because  little  effort  is 
required  to  split  the  problem  into  parallel  tasks.  Several  high- 
energy  Monte  Carlo  packages  have  been  ported  onto  computer 
clusters  using  tools  such  as  MPI  (Refs.  7,  8,  and  9)  and  shell 
scripting.10  Light  transport  in  turbid  medium  has  also  been  ac- 
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celerated  using  computer  clusters,11  multiprocessor  systems,12 
and  field-programmable  gate  arrays.13 

In  recent  years,  the  graphics  processing  unit  has  be¬ 
come  a  popular  platform  for  running  distributed  biomedical 
computations.14  For  simple  homogeneous  media,  GPU  com¬ 
puting  can  dramatically  accelerate  Monte  Carlo  simulation  of 
photon  migration.15  Very  high  acceleration  can  also  be  achieved 
for  voxelized  and  multilayer  geometries.16,17  Acceleration  is 
more  modest  for  complex  mesh-based  geometries18  because 
the  distributed  calculation  of  the  intersection  of  a  set  of  rays 
with  a  triangular  mesh  is  challenging  on  the  GPU.  Monte  Carlo 
simulation  for  high-energy  physics  has  also  been  investigated 
on  the  GPU.19,20  Because  high-energy  particles  can  undergo 
a  wide  range  of  physical  interactions  with  matter,  these  im¬ 
plementations  use  complex  optimization  strategies  for  efficient 
processing.21 

While  parallel  processing  techniques  can  accelerate  Monte 
Carlo  simulations,  practical  considerations  can  be  an  obstacle 
to  porting  existing  single-threaded  codes  onto  parallel  architec¬ 
tures.  To  utilize  parallel  resources  efficiently,  programmers  must 
be  skilled  in  parallel  programming  and  spend  substantial  effort 
optimizing  the  parallel  portion  of  their  code.  Parallel  code  is 
also  harder  to  debug  and  maintain.  Last,  large  computer  clus¬ 
ters  are  not  always  available  at  a  medical  institution  for  running 
massively  parallel  applications.  All  these  practical  issues  are 
important  drawbacks  to  the  development  and  use  of  distributed 
algorithms. 

Recently,  Internet- scale  computation  has  emerged  as  a  major 
driver  for  new  parallel  processing  technologies.  Internet  com¬ 
panies  routinely  process  very  large  datasets  such  as  log  files, 
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user  information,  pictures,  or  videos.  Most  of  these  processing 
tasks  are  quite  simple;  however,  the  size  of  the  input  data  is 
large  and  the  computation  has  to  be  distributed  to  hundreds  of 
nodes  for  practical  processing  times.  Usually,  porting  these  sim¬ 
ple  tasks  onto  parallel  architectures  requires  far  more  effort  than 
required  for  a  single-thread  implementation.  To  address  this  is¬ 
sue,  MapReduce  was  developed  at  Google  as  a  new  framework 
to  facilitate  the  development  of  parallel  algorithms.22 

MapReduce  can  hide  the  complexity  of  parallelization,  data 
distribution,  fault- tolerance,  and  load  balancing  to  the  developer, 
which  can  focus  on  developing  the  actual  algorithm.  In  this  pro¬ 
gramming  model,  the  developer  specifies  simple  tasks  which  are 
applied  in  a  distributed  fashion  to  large  datasets.  MapReduce  is 
also  well  integrated  with  existing  commercial  cloud  comput¬ 
ing  infrastructure.  Cloud  computing  refers  to  the  outsourcing  of 
one’s  compute  resources  to  third-party  companies.  Cloud  com¬ 
puting  providers  offer  services  such  as  web-based  software,  data 
warehousing,  and  scalable  clusters  of  virtual  nodes.  In  a  cloud 
computing  environment,  hardware  resources  are  often  virtual: 
computer  nodes  can  be  allocated  on  demand  with  custom  spec¬ 
ifications,  and  can  migrate  from  one  physical  host  to  another.  In 
a  cluster,  the  number  of  nodes  can  also  be  scaled  in  real-time, 
according  to  demand  for  computation. 

MapReduce  and  cloud  computing  technologies  are  al¬ 
ready  widely  used  in  applications  such  as  web  crawling 
and  analytics,23,24  data  mining,25  machine  learning,26  and 
bioinformatics.27-30  In  this  paper,  we  investigate  the  use  of 
MapReduce  for  biomedical  Monte  Carlo  simulation.  As  a  proof- 
of-concept,  we  ported  an  existing  photon  migration  Monte  Carlo 
code  to  Hadoop,  an  open-source  MapReduce  framework,  and 
characterized  the  performance  of  this  new  implementation  in  a 
cloud  computing  environment. 

2  MapReduce 

MapReduce  is  a  programming  model  for  processing  large  data 
sets.22  In  this  framework,  the  user  specifies  two  functions  called 
Map  and  Reduce,  respectively.  Although  neither  function  is  ex¬ 
plicitly  parallel,  many  instances  of  these  functions  are  executed 
concurrently  by  the  framework.  Input  data,  stored  on  a  dis¬ 
tributed  storage  system,  are  split  by  MapReduce  into  chunks 


and  distributed  to  parallel  Map  tasks  for  processing  (Fig.  1). 
No  communication  is  possible  between  Map  tasks;  however,  the 
outputs  of  Map  tasks  can  be  combined  by  Reduce  tasks. 

Data  communication  between  Map  and  Reduce  tasks  is  han¬ 
dled  using  key /value  pairs  (KVPs).  In  MapReduce,  keys  and 
values  can  be  stored  in  any  format,  provided  that  keys  can  be 
compared  to  one  another  and  sorted.  While  processing  input 
data,  Map  tasks  emit  a  sequence  of  intermediary  records  format¬ 
ted  as  KVPs  (Fig.  1).  Intermediary  records  that  share  a  common 
key  are  assigned  to  the  same  Reduce  task.  A  partition  function 
determines  how  keys  are  assigned  to  Reduce  tasks. 

The  role  of  the  Reduce  step  is  to  combine  intermediary 
records  into  the  final  output.  Each  Reduce  task  sequentially 
reads  records  associated  with  a  given  key.  The  user  specifies 
how  these  records  are  combined.  For  example,  the  user  can  sum 
the  values  associated  with  a  common  key.  The  output  of  the 
Reduce  function  is  automatically  written  in  KVP  format  to  a 
distributed  storage  system,  from  which  it  can  be  retrieved. 

Conceptually,  Map  and  Reduce  functions  can  be  described 
as 

map:  Vi  ^list(^2,v2) 

reduce:  (^2,list(v2))  ->  list(v3).  (1) 

The  number  of  Map  and  Reduce  tasks  is  set  by  the  user.  In  gen¬ 
eral,  the  number  of  tasks  should  be  larger  than  the  number  of 
nodes  for  fine  processing  granularity.  The  user  can  also  spec¬ 
ify  how  input  data  are  split  and  how  intermediary  records  are 
partitioned. 

A  feature  of  MapReduce  is  its  tolerance  for  hardware  fail¬ 
ures.  In  a  large  cluster,  the  probability  that  one  or  more  nodes 
fail  is  non-negligible.  When  a  node  fails,  MapReduce  resched¬ 
ules  the  compromised  tasks  onto  other  nodes  in  the  cluster.  A 
related  issue  is  the  heterogeneous  performance  of  worker  nodes. 
For  instance,  a  node  in  the  cluster  might  be  crippled  by  a  defec¬ 
tive  hardware  component  and  run  slowly.  Such  straggling  nodes 
can  substantially  reduce  the  performance  of  the  entire  cluster. 
MapReduce  uses  a  technique  called  speculative  execution  to 
overcome  this  type  of  situation.  Worker  nodes  that  finish  their 
workload  early  can  attempt  to  execute  clones  of  tasks  in  progress 


Write 


output  0 


>  output  1 


Output  files 
(HDFS) 


Fig.  1  An  overview  of  the  MapReduce  framework.  A  master  node  assigns  Map  and  Reduce  tasks  to  worker  nodes.  Input  files  are  split  into  chunks, 
that  are  processed  by  independent  Map  tasks,  producing  a  stream  of  intermediary  key/value  records.  These  records  are  selectively  read  by  Reduce 
tasks  according  to  their  key,  and  combined  finally  into  multiple  outputs. 
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Fig.  2  A  depiction  of  Hadoop  Streaming,  a  utility  that  allows  users  to 
specify  external  applications  as  Map  and  Reduce  functions. 

elsewhere  in  the  cluster.  A  task  is  considered  completed  when 
either  its  original  or  clone  completes. 

The  original  MapReduce  software  framework  developed  by 
Google  is  not  publicly  available  but  several  open-source  alter¬ 
natives  exist.  Hadoop,31  maintained  by  the  Apache  Software 
Foundation,  is  used  by  Internet  companies  for  large-scale  data 
processing,  financial  simulations  and  bioinformatics  calcula¬ 
tions.  Hadoop  set  a  record  in  2009  by  sorting  1  Petabyte  of  data 
in  16.3  h  on  3658  nodes.32 

The  Hadoop  project  includes  the  Hadoop  distributed  file  sys¬ 
tem  (HDFS),  designed  for  storing  extremely  large  data  files 
(Petabytes  and  up)  on  a  distributed  network  of  computers,  and 
Hadoop  MapReduce,  the  parallel  computation  engine.  Although 
Hadoop  is  written  in  Java,  developers  can  write  jobs  in  any  other 
programming  language  using  a  utility  called  Hadoop  Streaming. 
Hadoop  Streaming  implements  Map  and  Reduce  functions  as 
interfaces  to  external  user-specified  applications  (Fig.  2).  Exter¬ 
nal  Map  and  Reduce  applications  communicate  with  Hadoop 
Streaming  through  standard  Unix  streams.  They  read  input 
KVPs  via  standard  input  (stdin)  and  write  back  their  output 
via  standard  output  (stdout).  KVPs  can  be  formatted  as  Text  or 
TypedBytes,  which  are  sequences  of  bytes  in  which  the  first  byte 
is  a  type  code. 


3  Methods 

3.1  Monte  Carlo 

MC321,  a  steady-state  optical  Monte  Carlo  package  for  simu¬ 
lating  photon  migration  in  homogeneous  turbid  media,2  was  im¬ 
plemented  using  Hadoop.  In  this  package,  photons  are  launched 
in  packets  with  a  virtual  weight  initially  set  to  1 .  As  the  photon 
packet  propagates  in  the  medium,  it  can  scatter  or  be  absorbed. 
The  step  size  is  sampled  randomly  based  on  the  probability  den¬ 
sity  function  of  free  path.  Once  the  photon  packet  has  reached 
an  absorption  site,  a  fraction  of  its  weight  is  absorbed  and  scored 
into  a  256x256  detection  grid.  Next,  the  photon  packet  direc¬ 
tion  is  updated  by  sampling  the  deflection  angle  in  a  manner 
consistent  with  anisotropic  light  scattering.  Once  the  weight  of 
a  photon  packet  drops  below  0.01,  a  Russian  roulette  test  with 
a  10%  survival  probability  is  performed.  If  the  packet  survives, 
its  weight  is  renormalized  and  the  simulation  continues.  If  the 
roulette  test  fails,  the  packet  is  terminated  and  a  new  photon 
packet  may  be  initiated.  No  index  mismatch  or  boundary  con¬ 
ditions  are  considered  in  this  simple  simulation.  Photons  that 
reach  the  phantom’s  boundary  escape  the  medium  and  are  auto¬ 
matically  terminated. 

The  MC321  package  is  a  simplified  Monte  Carlo  method 
that  has  been  utilized  to  study  treatment  planning  for  photody¬ 
namic  therapy2  and  x-ray  luminescence  imaging.33  MC321  can¬ 
not  model  heterogeneous  media  with  spatially-varying  optical 
properties,  or  photon  propagation  through  material  boundaries. 
However,  MC321  shares  common  features  with  more  accurate 
packages,  such  as  MCML  (Ref.  3)  or  tMCimg.5 

3.2  MapReduce  Implementation 

We  propose  two  novel  distributed  implementations  of  MC321 
using  MapReduce.  The  first  method,  termed  Monte  Carlo 
event-based  processing  (MC-EP),  splits  the  production  and 
aggregation  of  photon  histories  into  Map  and  Reduce  steps, 
respectively.  In  this  scheme,  Map  tasks  run  parallel  Monte 


map 

(extern.) 


reduce 

(extern.) 


Fig.  3  An  overview  of  the  MC-EP  implementation.  (Map)  For  each  random  seed  loaded  from  stdin,  N  photon  packets  are  simulated.  Individual 
photon  events  are  written  to  stdout.  (Reduce)  The  stream  of  photon  events,  now  sorted  by  key,  are  aggregated  such  that  the  total  photon  weight 
absorbed  at  every  location  is  computed  and  written  out  to  stdout. 
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Carlo  simulations,  emitting  intermediary  KVPs  every  time  a 
photon  packet  deposits  a  fraction  of  its  weight  into  the  medium 
(Fig.  3,  top  panel).  In  these  intermediary  records,  the  key  is  an 
integer  index  to  the  nearest  two  dimensional  (2D)  detection  bin, 
and  the  value  a  floating-point  number  representing  the  photon 
weight  absorbed  in  that  bin.  These  intermediary  records  are 
then  aggregated  by  parallel  Reduce  tasks  in  the  following  way: 
records  that  share  the  same  key  are  sequentially  read  while 
the  corresponding  values  are  accumulated  in  a  local  variable. 
Once  all  the  records  for  a  given  key  have  been  processed,  the 
Reduce  tasks  outputs  the  cumulated  photon  weight  absorbed  in 
the  corresponding  detection  bin  (Fig.  3,  bottom  panel). 

The  MC-EP  approach  is  highly  flexible  because  it  separates 
the  production  and  the  aggregation  of  photon  histories.  In  gen¬ 
eral,  little  modification  is  required  to  implement  the  MC-EP 
scheme  using  an  existing  single- threaded  Monte  Carlo  code. 
Streaming  photon  events  to  stdout  using  the  KVP  format  re¬ 
quires  minimal  code  change.  (Typically,  only  a  few  fwrite  state¬ 
ments  need  to  be  added  to  the  code.)  However,  a  drawback  of  the 
MC-EP  approach  is  that  large  amounts  of  intermediary  data  are 
produced.  To  decrease  the  amount  of  data  exchanged  between 
Map  and  Reduce  steps,  intermediary  KVPs  can  be  partially  ag¬ 
gregated  within  the  Map  stage  prior  to  being  transferred  to  the 
Reduce  stage. 

In  this  scheme,  termed  Monte  Carlo  histogram-based  pro¬ 
cessing  (MC-HP),  each  Map  task  outputs  a  partial  scoring  array 
corresponding  to  the  photon  histories  it  has  computed.  In  turn, 
Reduce  tasks  accumulate  these  partial  outputs  into  a  complete 
scoring  array.  In  the  MC-HP  scheme,  intermediary  KVPs  consist 
of  a  unique  key  (arbitrarily  set  to  1)  and  an  array  of  floating-point 
values  representing  partial  2D  scoring  arrays.  Most  Monte  Carlo 
packages  (including  MC321)  score  photon  events  using  a  local 
array,  which  they  write  to  a  disk  once  the  simulation  completes. 
To  implement  the  MC-HP  scheme,  the  developer  only  needs  to 
write  the  scoring  array  to  stdout  instead  of  a  local  file. 

For  both  implementations,  Mapper  and  Reducer  are  written 
as  standalone  C  applications.  To  ensure  that  Monte  Carlo  simu¬ 
lations  performed  in  parallel  are  uncorrelated,  a  text  file  contain¬ 
ing  10,000  distinct  seeds  is  used  as  the  input  of  the  MapReduce 
job.  The  large  number  of  seeds  allows  Hadoop  to  execute  up  to 
10,000  parallel  Map  tasks  for  fine  processing  granularity.  We 
also  ensure  that  the  Map  function,  given  the  same  input,  always 
produces  the  same  output.  Random  number  generators  (RNGs) 
that  use  a  system-wide  state  are  problematic  with  MapReduce 
because  re-execution  and  speculative  execution  assume  deter¬ 
ministic  functions.  All  communications  through  UNIX  streams 
use  the  TypedBytes  data  format. 

The  MapReduce  implementation  of  the  MC321  Monte 
Carlo  package  described  in  this  paper  is  available  for  research 
purpose.  The  codes  can  be  downloaded  from  http://xinglab. 
stanford.edu/research/downloads.html. 

3.3  Computing  Environment 

Two  different  Hadoop  clusters  were  used.  To  develop  and  de¬ 
bug  MapReduce  jobs,  we  set  up  a  pseudo-distributed  cluster 
by  installing  Hadoop  0.21  on  a  quad-core  computer.  Input  and 
output  data  were  stored  on  a  local  installation  of  HDFS.  To  run 
and  profile  large  jobs,  we  allocated  a  remote  Hadoop  cluster 
on  Amazon’s  Elastic  Compute  Cloud  (EC2),  using  the  Elas¬ 


tic  MapReduce  (EMR)  service.  EMR  is  an  implementation  of 
Hadoop  0.20  tightly  integrated  with  other  Amazon  web  services, 
such  as  the  Simple  Storage  Service  (S3). 

Map  and  Reduce  applications  were  compiled  remotely  on 
an  EC2  node  using  GCC  version  4.3.2,  and  uploaded  onto  S3 
together  with  the  input  random  seeds.  EMR  jobs  were  submitted 
from  a  local  computer  using  a  freely-available  command-line 
tool. 

3.4  System  Evaluation 

For  benchmarking  purposes,  we  simulated  the  diffusion  of  a 
633  nm  laser  beam  within  a  turbid  5 -cm  radius  cylinder  filled 
with  a  5%  milk/95  %  water  mixture.34  The  phantom  properties 
were  set  as  follows:  the  absorption  and  scatter  coefficients  were 
/za  =  0.0025  cm-1  and  /zs  =  7.8  cm-1,  respectively;  the  index 
of  refraction  was  n  =  1.3;  and  the  anisotropy  coefficient  was 
g  =  0.9.  The  photons  were  initialized  at  a  single  point  on  the 
surface  of  the  phantom  and  directed  toward  its  center.  For  sim¬ 
plicity,  only  the  v  and  y  coordinates  of  each  absorbed  photon 
were  scored. 

In  a  first  experiment,  we  compared  the  original  MC321  code 
against  the  MC-EP  and  MC-HP  implementations.  The  bench¬ 
mark  consisted  in  simulating  100  million  photon  packets.  The 
original  Monte  Carlo  code  was  run  single-threaded  on  an  Intel 
Core2  Q9650  CPU  with  4  GB  of  memory.  The  MapReduce  im¬ 
plementations  were  run  on  EC2  using  20  high-memory  nodes 
(code-named  m2.xlarge).  The  m2.xlarge  node  configuration  is 
specially  adapted  for  tasks  with  large  memory  requirements, 
which  we  found  was  necessary  for  the  MC-EP  implementa¬ 
tion.  A  single  m2.xlarge  node  comes  with  17.1  GB  of  memory, 
two  64-bit  virtual  cores  with  3.25  EC2  compute  units  each,  and 
420  GB  of  local  storage.  The  output  data  were  downloaded  back 
to  a  local  computer  and  analyzed  with  matlab. 

In  a  second  experiment,  the  total  simulation  time  was 
recorded  for  a  variable  number  of  photon  packets,  ranging  from 
10  million  to  100  billion.  At  this  point,  we  focused  our  inves¬ 
tigation  on  the  more  efficient  MC-HP  implementation.  Twenty 
high-CPU  nodes  (codename  cl. medium)  were  used  in  these 
tests.  These  nodes  have  1.7  GB  of  memory,  two  32-bit  virtual 
cores  with  2.5  EC2  compute  units  each,  and  350  GB  of  local 
storage.  Each  cl. medium  node  is  configured  for  running  four 
MapReduce  tasks  simultaneously. 

In  a  third  experiment,  100  billion  photon  packets  were  sim¬ 
ulated  using  MC-HP  on  a  variable  number  of  nodes,  ranging 
from  1  to  240.  The  cl. medium  node  configuration  was  used. 
Five  thousand  Map  tasks  and  one  Reduce  task  were  run,  provid¬ 
ing  very  fine  task  granularity  and  good  load  balancing.  The  total 
run  time  was  recorded  and  the  output  of  the  simulations  com¬ 
pared.  Because  we  were  not  authorized  to  allocate  more  than 
240  EC2  nodes,  we  additionally  ran  the  same  simulation  on 
the  high-CPU  extra-large  node  configuration  (cl.xlarge).  These 
nodes  have  8  virtual  cores  with  2.5  EC2  compute  unit  each, 
7  Gb  of  memory,  and  1.7  Tb  of  local  storage.  A  cluster  of  240 
cl.xlarge  nodes  is  roughly  equivalent  to  960  cl. medium  nodes. 

In  a  fourth  experiment,  a  4-node  Hadoop  cluster  was  allo¬ 
cated  on  EC2  and  profiled  using  built-in  tools  while  simulating 
300  million  photons.  The  number  of  Map  and  Reduce  tasks  was 
set  to  100  and  1,  respectively.  Disk  I/O,  network  I/O,  and  CPU 
utilization  were  recorded  every  minute. 
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Fig.  4  MapReduce  execution  timeline  for  the  MC-EP  and  MC-HP  implementations,  showing  the  number  of  active  tasks  at  each  time  point.  The 
MC-EP  timeline  spans  1 8  min  09  s,  versus  1  min  54  s  for  the  MC-HP  timeline.  Waste  denotes  speculative  execution. 


In  a  last  experiment,  100  million  photon  packets  were  simu¬ 
lated  with  MC-HP  on  20  nodes  of  type  cl. medium.  After  3  min, 
a  set  number  of  worker  nodes  were  terminated  to  mimic  node 
failure  and  test  Hadoop’s  fault  tolerance. 

4  Results 

The  simulation  time  for  100  million  photon  packets  is  re¬ 
ported  for  the  original  and  distributed  implementations  (Table  1). 
For  MapReduce,  processing  time  is  defined  as  the  duration  of 
the  Hadoop  job,  not  including  cluster  initialization.  For  both 
MapReduce  jobs,  a  timeline  was  generated  from  Hadoop  log 
files  using  Karmasphere  Studio  (Fig.  4).  Overall,  MC-EP  is 
slower  than  MC-HP  because  MC-EP  transfers  many  intermedi¬ 
ary  records  between  Map  and  Reduce  tasks.  In  MC-EP,  shuffle, 
merge,  and  reduce  tasks  represent  a  significant  portion  of  all 
tasks  (Fig.  4,  top).  In  contrast,  MC-HP  runs  a  single  reduce  task 
(Fig.  4,  bottom). 

The  output  of  the  MC-HP  simulation  is  compared  for  a  vari¬ 
able  number  of  simulated  photon  packets,  ranging  from  10  mil¬ 
lion  to  100  billion  [Figs.  5(a)-5(e)].  The  output  of  the  MC-HP 
simulation  is  virtually  identical  to  that  of  the  original  MC321 
package.  For  10  million  photon  packets,  the  relative  root-mean- 
square  error  between  the  two  photon  absorption  distributions  is 
lower  than  4.1  x  10  _  7  [Fig.  5(f)].  The  discrepancy  between  the 
two  outputs  is  due  to  round-off  errors  during  the  conversion  of 
the  Reduce  output  to  10  digit  ASCII  in  the  final  file.  The  sim¬ 
ulation  time  also  increases  linearly  with  the  number  of  photons 
packets  [Fig.  6(a)].  Note  that  the  curvature  of  the  linear  fit  is  due 
to  the  logarithmic  scale  on  both  axes. 


Table  1  Simulation  time  comparison. 


No.  Nodes 

No.  Maps  No.  Red. 

Node  Type 

Sim.  Time 

MC321 

1 

N/A 

N/A 

Q9650 

28  min  39  s 

MC-EP 

20 

200 

200 

m2.xlarge 

1  8  min  09  s 

MC-HP 

20 

200 

1 

m2.xlarge 

1  min  54  s 

For  100  billion  photon  packets,  the  simulation  time  scaled 
inversely  with  the  number  of  nodes  [Fig.  6(b)].  On  a  single  dual¬ 
core  EC2  node,  simulating  100  billion  photons  took  11.2  days, 
a  1.8  x  speed-up  compared  with  the  single-threaded  MC321 
implementation.  The  same  simulation  took  only  22  min  on 
240  cl.xlarge  nodes,  a  1258  x  speed-up.  Nodes  based  on  the 
cl.xlarge  (8  cores)  configuration  ran  the  MC-HP  simulation  3.3 
times  faster  than  cl. medium  nodes  (2  cores). 

To  better  understand  the  use  of  resources,  a  4-node  Hadoop 
cluster  was  profiled  while  simulating  300  million  photons 
(Fig.  7).  The  total  runtime  for  this  experiment  was  20  min.  It 
can  be  observed  that  the  execution  of  the  simulation  was  clearly 
compute-bound  because  the  CPU  utilization  was  maximized  on 
all  three  slave  nodes.  Only  the  CPU  capacity  of  the  master  node 
was  not  fully  utilized.  The  pattern  of  network  I/O  was  different 
for  nodes  running  Map  tasks  only  and  nodes  running  both  Map 
and  Reduce  tasks.  In  MC-HP,  Map  tasks  are  data  producers;  they 
input  little  data  (the  random  seeds)  but  output  large  amounts  of 


(a)  1(b)  1(c) 


10M  100M  IB 


10B  100B  Error 


Fig.  5  (a)-(e)  Output  of  the  MC-HP  simulation  for  a  variable  number 
of  photon  histories  ranging  from  1 0  million  to  1 00  billion,  shown  on  a 
base-10  logarithmic  scale,  (f)  Error  between  the  original  Monte  Carlo 
package  and  the  MapReduce  implementation,  shown  for  10  million 
photon  histories  on  a  base-10  logarithmic  scale. 
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Fig.  6  (a)  MC-HP  simulation  time  for  a  varying  number  of  photon 
histories.  Dashed  line:  least-square  linear  regression  ( r 2  =  0.9998).  The 
curvature  of  the  linear  fit  is  due  to  the  logarithmic  scale,  (b)  MC-HP 
simulation  time  for  a  varying  number  of  CPU  cores.  Dashed  line:  ideal 
inverse  scaling  of  the  single-node  simulation  time  (r2  =  0.98). 


data  (the  scoring  arrays).  In  contrast,  Reduce  tasks  are  data  con¬ 
sumers:  they  continuously  receive  data  from  Map  tasks,  which 
they  only  write  to  HDFS  at  the  end  of  the  simulation. 

Hadoop’ s  fault  tolerance  was  also  evaluated  by  shutting  down 
1  and  10  worker  nodes,  respectively,  during  the  simulation  of 
100  million  photons  on  20  nodes.  The  termination  of  these  nodes 
did  not  prevent  Hadoop  from  successfully  completing  its  job, 
provided  that  the  master  node  was  kept  alive.  Upon  node  failure, 
EMR  allocated  new  nodes  dynamically  and  relaunched  failed 
Map  and  Reduce  tasks.  The  output  of  the  Hadoop  simulation 
with  some  nodes  terminated  was  identical  to  that  obtained  by 
running  the  simulation  with  no  failure.  The  failure  of  1  and 
10  nodes  lengthened  the  total  simulation  by  15  and  18  min, 
respectively. 


5  Discussion 

While  MapReduce  is  often  used  to  crunch  large  datasets  such 
as  website  logs,  we  demonstrated  that  it  is  also  well  suited  for 
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Fig.  7  Resource  profile  showing  I/O  and  computing  during  simulation 
of  300  million  photons  on  a  4-node  Hadoop  cluster. 


computationally-demanding  Monte  Carlo  simulations.  MapRe¬ 
duce’s  architecture  is  largely  decentralized  for  efficient  data  pro¬ 
cessing  on  large  clusters.  Even  with  240  nodes,  the  performance 
of  the  MC-HP  implementation  scales  linearly  with  the  number 
of  CPU  cores  [Fig.  6(b)],  showing  few  signs  of  diminishing  re¬ 
turns.  On  average,  the  simulation  times  were  within  12%  of  the 
ideal  inverse  scaling  [Fig.  6(b),  dashed  line].  This  result  suggests 
that,  if  needed,  the  MC-HP  implementation  could  run  on  many 
more  nodes  without  suffering  from  interconnect  latencies  and 
limited  network  throughput. 

The  MC-HP  simulation  also  scaled  linearly  with  the  num¬ 
ber  of  photon  packets  simulated.  When  the  number  of  photons 
was  incremented  from  10  million  to  100  billion,  the  simulation 
time  followed  a  linear  trend  [Fig.  6(a),  r2  =  0.9998].  Each  in¬ 
cremental  one  million  photons  can  be  simulated  on  20  nodes  in 
0.6  s.  Furthermore,  the  computation  overhead,  which  consists  of 
cluster  initialization,  task  distribution,  and  data  reduction,  is  on 
the  order  of  100  s,  which  is  negligible  when  billions  of  photon 
packets  are  simulated  [Fig.  6(a)]. 

By  combining  the  information  from  the  first  and  second  ex¬ 
periment,  the  performance  of  the  MC-HP  implementation  can 
be  better  characterized.  The  throughput  and  latency  were  found 
to  be  85,178  photons  per  node  per  second,  and  100  s,  respec¬ 
tively.  In  a  separate  experiment  (data  not  shown),  latency  was 
found  to  be  independent  of  the  number  of  nodes. 
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Hadoop  also  proved  to  be  a  highly  reliable  framework  that 
could  withstand  the  accidental  failure  of  50%  of  the  worker 
nodes.  Upon  node  failure,  EMR  allocates  new  nodes  to  replace 
failed  ones.  Currently,  the  process  of  resizing  a  running  Hadoop 
cluster  on  EC2  is  rather  slow,  requiring  around  15  min,  but  a 
substantial  effort  is  under  way  to  improve  the  efficiency  of  the 
process.  There  also  exist  simpler  implementations  of  MapRe¬ 
duce  that  do  not  resize  the  cluster,  but  rather  redistribute  failed 
tasks  to  nodes  not  affected  by  the  failure. 

A  drawback  of  cloud  computing  is  that  a  virtual  cluster  must 
be  allocated  on  the  cloud  before  running  a  job,  a  process  that  in¬ 
volves  deploying  a  large  number  of  operating  system  images  on 
physical  computers.  Such  initialization  contributes  to  the  over¬ 
all  execution  time  by  increasing  the  latency,  especially  for  very 
large  cluster  configurations.  The  solution  we  adopted  consisted 
in  keeping  our  virtual  Hadoop  cluster  alive  after  the  completion 
of  its  first  job,  such  that  it  would  remain  available  for  other  jobs. 

A  word  should  also  be  said  about  the  cost  of  performing 
Monte  Carlo  simulation  in  a  commercial  cloud  computing  envi¬ 
ronment.  As  of  August  2011,  the  price  of  Amazon’s  EC2  nodes 
ranges  from  $0.025/h  for  micronodes  to  $2.28/h  for  quadru¬ 
ple  extra-large  nodes.  Most  of  the  simulations  presented  in  our 
study  were  performed  on  high-CPU  medium  nodes,  which  are 
priced  at  $0.1 9/h.  Hence,  large-scale  simulations,  such  as  those 
with  1011  photon  packet  histories  [Fig.  6(b)],  can  cost  a  few 
hundred  dollars  to  run.  Cloud  computing  is  ideal  for  research 
institutions  that  do  not  have  a  sufficient  demand  to  maintain 
a  full-size  cluster  on  site.  Furthermore,  cloud  resources  have 
very  low  downtime.  The  main  alternatives  to  cloud  computing 
is  the  operation  of  a  dedicated  cluster,  a  solution  currently  used 
by  many  research  groups.  This  approach  has  a  few  advantages, 
namely  the  absence  of  a  virtualization  layer  and  better  control 
over  the  hardware.  However,  energy,  maintenance,  and  storage 
costs  can  be  significantly  higher.  Furthermore,  the  lack  of  elastic 
scaling  results  in  under-  or  overutilization  of  compute  resources. 

MapReduce  is  highly  optimized  for  data-intensive  process¬ 
ing  tasks,  such  as  indexing  web  pages  or  mining  large  datasets. 
It  uses  a  loosely-coupled  distributed  memory  system,  which 
does  not  allow  parallel  tasks  to  communicate  with  one  another 
directly.  By  its  design,  MapReduce  is  suboptimal  for  compute- 
bound  tasks  or  for  tasks  that  require  frequent  internode  data 
exchange,  such  as,  for  instance,  iterative  image  reconstruction. 
In  MapReduce,  internode  communication  can  be  achieved  by 
reducing  the  Map  outputs  with  Reduce  tasks,  and  splitting  these 
data  again  using  a  new  set  of  Map  tasks.  For  tasks  with  a  large 
amount  of  internode  communication,  this  mechanism  may  be 
less  efficient  than  a  shared  memory  approach.  For  instance, 
GPUs  provide  on-chip  memory  shared  by  threads  in  a  same 
block,  and  global  memory  shared  by  all  threads.  As  a  result, 
GPU  threads  can  exchange  data  efficiently.  Many  scientific  com¬ 
putation  strategies  have  used  the  GPU’s  fast  shared  memory  to 
achieve  impressive  performance.33  Furthermore,  GPU  can  ac¬ 
cess  global  memory  with  nearly  5 -fold  higher  bandwith  than 
state-of-the-art  CPUs:  For  instance,  the  memory  bandwidth  of 
a  Tesla  C2050  GPU  is  144  GB/s,  versus  32  GB/s  for  an  In¬ 
tel  Xeon  X5650  CPU.  However,  the  peak  performance  of  the 
two  architectures  should  be  compared  with  caution  because  the 
GPU’s  peak  memory  bandwidth  is  only  achieved  when  mem¬ 
ory  transfers  are  fully  coalesced  (i.e.,  they  occur  in  continous 
128-bit  chunks).  For  random  memory  accesses  (such  as  scoring 
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during  a  Monte  Carlo  simulation),  this  is  seldom  the  case  and  the 
achieved  bandwith  can  be  lower  than  the  peak  theoretical  value. 
At  this  time,  the  use  of  a  distributed  memory  system  such  as 
MapReduce  for  scientific  applications  other  than  Monte  Carlo 
is  still  being  investigated.35 

The  performance  of  the  MC-HP  implementation  can  be  put 
into  perspective  by  comparing  it  against  other  published  simu¬ 
lations  of  photon  migration  in  homogeneous  media  (Table  2). 
These  performance  figures  should  be  compared  with  caution 
because  of  the  diverse  set  of  parameters  used  in  these  studies. 
The  speed  of  the  simulation  can  be  affected  by  a  number  of  pa¬ 
rameters,  including  the  RNG,  the  optical  properties,  the  size  of 
the  scoring  grid,  and  the  roulette  parameters.  With  this  caveat  in 
mind,  the  cloud-based  approach  provides  speed-up  of  the  same 
order  of  magnitude  as  its  GPU-based  counterparts  (Table  2).  A 
GPU  platform  has  several  advantages  for  Monte  Carlo  simula¬ 
tion,  such  as  built-in  hardware-accelerated  arithmetic  functions, 
no  thread-switching  overhead,  a  compact  form  factor,  and  the 
relatively  low  price  of  the  hardware.14 

However,  one  issue  with  GPU  computing  is  that  many  threads 
share  access  to  the  same  memory.  Hence,  race  conditions  may 
arise  when  two  or  more  threads  attempt  to  write  to  the  same 
grid  element,  since  only  one  thread  can  write  to  a  given  mem¬ 
ory  location  on  a  clock  cycle.  This  situation  occurs  most  fre¬ 
quently  in  high  fluence  regions,  such  as  near  the  light  source, 
and  is  less  problematic  in  low-fluence  regions.16  GPU  atomic 
operations  can  be  used  to  serialize  the  memory  writes  when 
needed,  at  the  cost  of  slower  memory  access.  For  instance,  the 
MCX  package  is  four  times  slower  when  such  atomic  opera¬ 
tions  are  used  (Table  2).  A  variant  of  MCX  (mcx_cached)  and 
the  GPU-MCMF  package  use  a  different  approach  to  mitigate 
race  condition  issues:  The  high-fluence  region  of  the  scoring 
array  is  cached  in  fast  shared  memory  and  accessed  atomically, 
while  the  low-fluence  region  is  stored  in  global  memory.  Use  of 
fast  shared  memory  reduces  the  penalty  associated  with  atomic 
operations  since  fewer  threads  access  the  grid  at  any  time,  and 
shared  memory  has  much  lower  latency  than  global  memory. 
In  our  MapReduce  implementation,  each  task  allocates  its  own 
detection  grid  locally,  which  avoids  data  write  hazards.  Focal 
detection  grids  are  combined  using  parallel  Reduce  tasks. 

Additionally,  MapReduce  offers  better  scalability  than 
GPUs.  To  achieve  higher  photon  throughput,  multiple  GPUs 
must  be  employed  and  integrated  into  a  cluster.17  The  software 
must  be  adapted  to  provide  inter- GPU  and  intercomputer  com¬ 
munications,  which  adds  another  layer  of  complexity.  In  con¬ 
trast,  the  number  of  nodes  in  a  MapReduce  cluster  can  be  scaled 
by  modifying  a  single  parameter.  However,  if  used  intensively, 
a  GPU  cluster  can  be  more  economical  over  time  than  com¬ 
puting  capacity  purchased  on  demand  from  a  cloud  computing 
operator. 

In  the  context  of  a  Monte  Carlo  simulation,  a  loosely-coupled 
parallel  architecture  such  as  a  Hadoop  cluster  has  some  ad¬ 
vantages  over  a  GPU.  The  GPU  architecture  is  based  on  a 
single-program  multiple- data  programming  model  which  pe¬ 
nalizes  diverging  threads.  Briefly,  a  GPU  is  composed  of  several 
multiprocessors,  each  of  which  processes  parallel  threads  in 
batches  called  warps.  Even  though  threads  within  a  warp  can 
follow  diverging  paths,  a  GPU  multiprocessor  can  only  issue 
one  instruction  per  clock  cycle.  Hence,  the  execution  of  di¬ 
verging  threads  is  serialized:  only  those  threads  programmed  to 
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Table  2  Comparison  of  various  photon  migration  Monte  Carlo  codes,  according  to  the  following  parameters:  RNG;  the  absorption  and  scatter 
coefficients  /za  and  /zs,  respectively;  the  size  of  the  scoring  grid;  the  acceleration  platform;  the  photon  throughput  Tacc  and  Tref,  expressed  in 
photons/ms,  for  the  accelerated  and  reference  (CPU)  platform,  respectively;  and  the  speed  up  achieved. 


Implementation 

RNG 

/xatcrrn1) 

^s(cm_1 ) 

Scoring  grid 

Platform 

Threads 

face 

^"ref 

Speed-up 

MC-HP 

RAN  3 

0.0025 

7.8 

256  x  256 

EC2  cloud 

3,840 

75,471 

60 

1258 

MCX  (atomic)  (Ref.  1 6) 

LL5 

0.05 

10.1 

O' 

o 

X 

O' 

o 

X 

O' 

o 

G92 

512 

900° 

12 

75 

MCX  (approx.)  (Ref.  16) 

LL5 

0.05 

10.1 

60  x  60  x  60 

G92 

1792 

3,800° 

12 

325 

WMC-GPU  (Ref.  15) 

MWC 

0 

90 

200a 

8800GT 

26,880 

1,747 

1.6 

1080 

GPU-MCML  (Ref.  1 7) 

MWC 

0.015 

708 

1 00  x  1 00 

GTX  480 

13,440 

384 

0.4 

870 

a  Estimated. 


execute  the  current  instruction  can  advance  to  the  next  in¬ 
struction;  others  must  wait  for  the  multiprocessor  to  execute 
their  instruction.  In  a  Monte  Carlo  simulation,  particle  paths 
can  diverge.  For  instance,  high-energy  photons  can  undergo  a 
wide  range  of  physical  interactions  with  matter.  Furthermore, 
a  comprehensive  model  of  dose  deposition  in  matter  requires 
simulation  of  the  electrons  induced  by  Compton  scatter  and 
photoelectric  absorption.  Because  of  these  effects,  naive  GPU 
implementations  result  in  high  thread  divergence  and  modest 
efficiency.10,20  Simulation  of  photon  migration  in  turbid  media 
can  also  lead  to  thread  divergence:  within  one  iteration,  some 
photons  may  traverse  a  boundary  or  be  reflected;  others  may  be 
scattered  or  absorbed.  Unlike  GPU  cores,  CPU  cores  have  dedi¬ 
cated  instruction  units.  Therefore,  each  parallel  thread  can  issue 
its  own  instruction  on  a  clock  cycle,  independently  of  the  other 
threads.  In  a  MapReduce  job,  parallel  tasks  are  executed  as  in¬ 
dependent  processes  by  CPU  cores,  and  are  not  penalized  if  they 
diverge  and  follow  different  code  branches.  It  should  be  noted 
that  the  cores  in  a  modern  CPU  can  also  issue  special  streaming 
SIMD  extension  (SSE)  instructions,  which  are  applied  simul¬ 
taneously  to  up  to  four  single-precision  floating-point  registers. 
This  lower  level  of  parallelism  cannot  be  automatically  utilized 
by  MapReduce  and  requires  the  developer  to  specifically  use 
SSE  instructions  when  developping  the  Map  and  Reduce  tasks. 

The  Monte  Carlo  code  used  in  this  paper  was  intended  as 
a  case  study  of  MapReduce  for  a  biomedical  application.  The 
simplicity  of  the  Monte  Carlo  code  helped  focus  our  presen¬ 
tation  on  the  mechanisms  of  task  distribution.  The  principles 
described  in  this  paper  can  be  applied  to  more  complex  photon 
migration  simulations,  such  as  those  that  include  a  three  dimen¬ 
sional  (3D)  map  of  the  optical  coefficients5  or  a  multilayered 
geometry.3  Hadoop  has  a  distributed  cache  that  can  make  a  set 
of  files  available  to  all  the  nodes  in  the  cluster.  This  mecha¬ 
nism  can  be  used  to  broadcast  the  optical  coefficient  map  to 
all  the  nodes  prior  to  starting  the  simulation.  In  our  implemen¬ 
tation,  the  distributed  cache  was  used  to  transfer  the  MC321 
binaries  (100  kB)  to  all  the  nodes.  Once  the  optical  coefficients 
are  available  on  all  the  nodes,  complex  Monte  Carlo  simulations 
using  heterogeneous  geometries  may  be  performed  by  parallel 
map  tasks,  with  no  internode  communications.  The  outputs  of 
these  parallel  simulations  may  be  combined  by  simple  reduced 
tasks,  as  demonstrated  in  this  work.  Apart  from  the  initial  broad¬ 
cast  of  the  optical  coefficient  map,  the  use  of  a  more  accurate 


Monte  Carlo  package  would  require  no  additional  data  transfer 
between  nodes.  The  present  work  suggests  that  these  more  com¬ 
plex  Monte  Carlo  simulations  would  be  accelerated  by  a  factor 
approximately  equal  to  the  number  of  nodes  in  the  cluster. 

In  the  MC-HP  implementation,  a  singe  Reduce  task  is  used 
to  aggregate  the  outputs  from  all  Map  tasks.  This  approach  may 
not  be  fast  enough  to  reduce  large  outputs  (such  as  for  instance 
a  3D  photon  fluence  matrix)  from  many  Map  tasks.  Network 
throughput  may  also  be  problematic  since  data  are  transferred 
from  many  nodes  to  a  single  node,  creating  a  bottleneck.  To 
alleviate  this  issue,  parallel  Reduce  tasks  can  be  generated  by 
setting  Map  tasks  to  output  their  data  block  by  block,  using  the 
key  to  encode  the  block  index.  The  blocks  may  be  composed  of 
slices  for  a  3D  matrix,  or  rows  in  a  2D  matrix.  This  data  transfer 
strategy  has  the  advantage  that  it  is  decentralized:  many  data 
transfer  can  occur  in  parallel,  increasing  the  total  bandwidth  of 
the  cluster.  For  optimal  efficiency,  the  number  of  Reduce  tasks 
can  be  set  to  the  number  of  such  blocks.  After  the  Reduce  step, 
the  blocks — which  combine  the  computations  of  many  Map 
tasks — are  written  to  distributed  storage,  from  which  they  can 
be  downloaded. 

Although  it  has  its  roots  in  text  processing,  Hadoop  is  an 
evolving  technology,  continuously  adding  new  features  to  ad¬ 
dress  the  specific  needs  of  a  wider  spectrum  of  users.  As  the 
framework  evolves,  we  expect  that  it  will  become  more  general- 
purpose,  and  even  better  suited  for  scientific  computation.  For 
example,  Hadoop  currently  requires  worker  nodes  to  sort  inter¬ 
mediary  records  by  key  before  applying  the  Reduce  function. 
This  sorting  step  is  unnecessary  in  scientific  applications  that 
use  keys  as  array  indices.  The  open-source  status  of  the  Hadoop 
project  allows  any  developer  to  modify  the  framework  as  needed, 
and  possibly  contribute  back  to  the  project. 

6  Conclusion 

The  inexorable  shift  of  computing  to  parallel  architectures  forces 
us  to  rethink  how  algorithms  are  implemented.  Porting  an  ap¬ 
plication  to  a  distributed  environment  is  problematic  because  of 
the  complexity  of  the  software  development  and  the  cost  of  a 
large-scale  computer  cluster.  For  photon  migration  Monte  Carlo 
simulation,  MapReduce  helps  overcome  these  two  challenges 
by  providing  a  simple  way  to  deploy  massively -parallel  applica¬ 
tions  in  a  cloud  computing  environment.  The  port  of  the  MC321 
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package  to  MapReduce  was  rapid,  taking  only  a  week;  yet,  the 
code  execution  scaled  to  1920  cores  with  few  signs  of  diminish¬ 
ing  returns.  The  new  distributed  code  includes  features  such  as 
fault  tolerance  and  automated  load  balancing. 
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Radioluminescence  microscopy  is  a  new  method  for  imaging  ra¬ 
dionuclide  uptake  by  single  live  cells  with  a  fluorescence  micro¬ 
scope.  Here,  we  report  a  particle-counting  scheme  that  improves 
spatial  resolution  by  overcoming  the  (3-range  limit.  Methods:  Short 
frames  (10  pus— 1  s)  were  acquired  using  a  high-gain  camera  cou¬ 
pled  to  a  microscope  to  capture  individual  ionization  tracks.  Optical 
reconstruction  of  the  (3-ionization  track  (ORBIT)  was  performed  to 
localize  individual  (3  decays,  which  were  aggregated  into  a  compos¬ 
ite  image.  The  new  approach  was  evaluated  by  imaging  the  uptake 
of  18F-FDG  in  nonconfluent  breast  cancer  cells.  Results:  After  im¬ 
age  reconstruction,  ORBIT  resulted  in  better  definition  of  individual 
cells.  This  effect  was  particularly  noticeable  in  small  clusters  (2-4 
cells),  which  occur  naturally  even  for  nonconfluent  cell  cultures.  The 
annihilation  and  Bremsstrahlung  photon  background  signal  was 
markedly  lower.  Single-cell  measurements  of  18F-FDG  uptake  that 
were  computed  from  ORBIT  images  more  closely  matched  the  up¬ 
take  of  the  fluorescent  glucose  analog  (Pearson  correlation  coef¬ 
ficient,  0.54  vs.  0.44,  respectively).  Conclusion:  ORBIT  can  image 
the  uptake  of  a  radiotracer  in  living  cells  with  spatial  resolution  better 
than  the  (3  range.  In  principle,  ORBIT  may  also  allow  for  greater  quan¬ 
titative  accuracy  because  the  decay  rate  is  measured  more  directly, 
with  no  dependency  on  the  (3-particle  energy. 

Key  Words:  radionuclide  imaging  instrumentation;  single-cell  anal¬ 
ysis;  microscopy;  autoradiography 
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A 

X  Autoradiography  is  a  well-established  method  for  high-resolution 
imaging  of  radionuclide  probes  in  tissues.  Film  and  emulsion  meth¬ 
ods  have  the  highest  spatial  resolution  but  poor  sensitivity,  dy¬ 
namic  range,  and  quantitative  accuracy  and  require  tedious  sample 
preparation  (1,2).  Other  autoradiography  methods  (e.g.,  storage 
phosphor  (3),  solid-state  detection  (4,5),  gaseous  detectors  (6),  and 
thin  phosphor  (7))  have  higher  sensitivity  and  dynamic  range  but 
spatial  resolution  worse  than  50  pan.  Only  a  few  methods  have 
demonstrated  the  ability  to  visualize  the  uptake  of  radionuclide 
probes  with  single-cell  resolution.  One  of  these  methods  uses  a  13- 
sensitive  avalanche  photodiode  to  measure  radionuclide  uptake  in 
small  groups  of  cells,  cultured  in  16  microfluidic  chambers  (8).  An 
experiment  showed  that  the  avalanche  photodiode  could  measure 
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signal  from  a  single  cell  in  the  chamber.  Another  device,  the  Micro- 
Imager,  achieved  15-pan  spatial  resolution  for  35S,  which  was  used 
to  detect  in  situ  hybridization  in  single  neurons  (9).  A  third  device, 
the  radioluminescence  microscope,  was  developed  to  visualize 
radionuclide  uptake  in  live  cells  during  fluorescence  microscopy 
(10).  Radioluminescence  microscopy  can  be  used  to  measure  fluo¬ 
rescent  and  radionuclide  signals  emanating  from  a  collection  of 
living  cells,  in  a  relatively  short  time.  Here  we  propose  a  new  particle¬ 
counting  scheme  for  radioluminescence  microscopy  with  higher 
spatial  resolution  and,  in  principle,  quantitative  accuracy. 

The  standard  radioluminescence  microscopy  approach  uses  a 
scintillator  plate,  placed  in  contact  with  the  cells  of  interest,  to 
convert  ionizing  radiation  from  radioactive  (3  decay  into  visible- 
range  photons  detectable  with  a  fluorescence  microscope  (Fig.  1A). 
To  ensure  good  contact,  the  cells  are  grown  such  that  they  adhere 
to  the  scintillator  plate.  The  transparency  of  the  plate  in  the  visible 
range  allows  for  conventional  microscopy  techniques  (e.g.,  fluo¬ 
rescence  and  brightfield).  Before  imaging,  the  cells  are  incubated 
with  a  radionuclide  probe  and  then  washed  to  remove  the  excess 
probe.  A  radioluminescence  image  is  acquired  using  a  sensitive 
camera  coupled  to  a  high-numerical- aperture  microscopy  objec¬ 
tive  (Fig.  1A).  A  typical  image  requires  an  exposure  time  of  sev¬ 
eral  minutes  and  combines  the  contributions  of  thousands  of  (3 
particles,  each  emitted  with  a  random  energy  and  direction. 

Although  radioluminescence  microscopy  can  measure  radio- 
tracer  uptake  with  single-cell  resolution  (10),  this  single-exposure 
approach  is  suboptimal  for  2  reasons.  First,  as  each  |3-particle  trav¬ 
els  through  the  scintillator  plate,  it  leaves  behind  a  track  of  optical 
photons.  Since  a  (3-particle  may  travel  beyond  the  boundaries  of 
the  cell  in  which  it  was  emitted,  it  may  contribute  radiolumines¬ 
cence  signal  to  other  surrounding  cells,  which  makes  it  difficult  to 
distinguish  individual  cells  in  the  radioluminescence  image.  Sec¬ 
ond,  the  contribution  of  each  individual  positron  decay  event  to 
the  total  luminescent  signal  varies  greatly  because  each  positron  is 
emitted  with  a  variable  amount  of  kinetic  energy  and  deposits  a 
variable  fraction  of  this  energy  in  the  scintillator  at  a  variable 
depth.  The  inclusion  in  the  signal  of  these  random  quantities  re¬ 
sults  in  higher  statistical  variance. 

We  report  a  new  image  acquisition  scheme  for  radioluminescence 
microscopy  that  solves  the  2  aforementioned  problems.  Rather  than 
acquiring  a  single  image  with  a  long  exposure  time  (e.g.,  5  min), 
we  take  thousands  of  short  frames  with  high-electron-multiplica¬ 
tion  gain.  With  this  method,  we  can  visualize  the  ionization  tracks 
produced  by  individual  (3  particles  as  they  travel  through  the  scin¬ 
tillator.  We  trace  each  track  backward  to  the  2-dimensional  loca¬ 
tion  of  the  corresponding  (3  decay.  From  the  estimated  locations  of 
many  radioactive  decays,  we  synthesize  a  composite  image  that 
has  improved  spatial  resolution  compared  with  the  single-exposure 
acquisition  scheme  (Fig.  IB).  We  have  given  this  particle-counting 
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FIGURE  1.  Radioluminescence  microscopy  using  ORBIT.  (A)  Emission  of  positron  from  radiotracer  bound  to  cell  (7)  is  followed  by  positron- 
electron  annihilation  (2)  and  emission  of  2  back-to-back  511-keV  photons  (blue  lines),  which  may  interact  with  electron  in  scintillator  via 
Compton  scatter  or  photoelectric  absorption  (3).  Energetic  charged  particles  (red  lines)  propagate  through  scintillator,  producing  light, 
which  can  be  recorded  by  microscope.  Because  of  shallow  depth  of  field,  system  is  sensitive  only  to  events  occurring  near  focal  plane. 
(Drawing  is  not  to  scale:  objective’s  field  of  view  is  1.2  mm).  (B)  To  acquire  radioluminescence  image  with  ORBIT,  bright-field  micrograph  is 
first  acquired  to  delineate  cell  boundaries.  Sequence  of  short  image  frames  is  then  acquired  with  high  gain  to  capture  ionization  tracks 
produced  by  energetic  charged  particles.  Those  frames  are  processed  to  extract  decay  locations,  which  are  then  aggregated  into  com¬ 
posite  image.  (Each  decay  location  is  represented  as  a  dot).  NA  =  numerical  aperture. 


scheme  the  name  optical  reconstruction  of  the  /3-ionization  track 

(ORBIT). 

Although  in  principle  similar  to  existing  particle-counting 
autoradiography  systems  (5,9),  ORBIT  can  image  and  reconstruct 
individual  ionization  tracks.  In  this  regard,  it  is  comparable  to 
single-molecule  fluorescence  microscopy  methods  such  as  subdif¬ 
fraction-limit  imaging  by  stochastic  optical  reconstruction  micros¬ 
copy  (11).  This  method  uses  bright  light  pulses  to  stochastically 
photoactivate  a  small  number  of  fluorescent  proteins,  which  are 
then  fluorescently  imaged.  Because  of  their  sparse  distribution, 
individual  activated  proteins  can  be  localized  with  precision  greater 
than  the  diffraction  limit.  By  repeating  this  process,  different  sub¬ 
sets  of  proteins  can  be  activated  and  imaged.  Similarly,  ORBIT 
exploits  the  fact  that  radioactive  decay  is  stochastic,  and  thus  only 
a  small  number  of  randomly  distributed  ionization  tracks  are  pro¬ 
duced  within  a  short  time  frame.  The  number  of  optical  photons 
produced  in  a  scintillator  during  a  single  radioactive  decay  event  is 
large  enough  to  be  detected  by  a  sensitive  microscope.  For  instance, 
the  kinetic  energy  of  a  positron  emitted  during  the  decay  of  18F 
can  be  as  high  as  634  keV,  which  for  a  scintillator  such  as  CdW04 
(15  photons/keV)  would  result  in  up  to  9,500  photons  emitted. 

MATERIALS  AND  METHODS 
Imaging  System 

The  microscope  setup  used  in  this  study  has  been  described  else¬ 
where  (10).  A  40x/1.3-numerical  aperture  oil-immersion  objective 
(UPLFLN40XO;  Olympus)  was  used  for  all  experiments.  We  estimate 
that  the  microscope  collects  approximately  13%  of  the  photons  emit¬ 
ted  near  the  focal  plane,  based  on  the  objective  light  collection  effi¬ 
ciency  (25%  for  a  1.3  numerical  aperture)  and  transmission  (80%),  the 
throughput  of  the  tube  lens  (80%),  and  the  quantum  efficiency  of  the 
charge-coupled  device  (CCD)  (80%).  The  light  collection  efficiency 
calculation  does  not  account  for  the  higher  index  of  refraction  of  the 
scintillator  (n  =  2.3).  Experimentally,  we  found  that  the  microscope 
could  record  31%  of  the  photons  emitted  by  a  bioluminescent  light 
source  (Supplemental  Fig.  1;  supplemental  materials  are  available 
online  at  http://jnm.snmjoumals.org).  The  discrepancy  between  these 


2  values  may  be  due  to  camera  calibration  factors,  electron-multipli¬ 
cation  gain  nonlinearity,  and  the  nonideal  behavior  of  the  microscope 
objective. 

Depth-Varying  Point-Spread  Function 

We  characterized  the  depth- varying  point-spread  function  and  mod¬ 
ulation  transfer  function  of  the  microscope  objective  using  an  analytic 
model  (12).  The  model  accounts  for  the  immersion  oil  and  the  index 
of  refraction  of  the  scintillator.  Using  this  model,  we  computed  via 
convolution  the  response  of  the  microscope  to  a  small  pointlike  object 
placed  at  different  depths.  The  object  was  a  3-dimensional  gaussian 
function  with  a  full  width  at  half  maximum  of  6.5  pm — that  is,  twice 
the  camera  pixel  size  (with  2x2  binning). 

ORBIT  Imaging 

The  general  methodology  for  acquiring  an  ORBIT  image  is  as  fol¬ 
lows.  We  first  acquire  a  sequence  of  frames  with  a  short  exposure 
time,  an  electron-multiplication  gain  of  1,200,  and  2x2  binning.  The 
exposure  time  must  be  set  to  maximize  the  system  sensitivity.  On  the 
one  hand,  the  exposure  should  be  short  enough  to  avoid  spatiotempo- 
ral  overlap  between  different  ionization  tracks.  On  the  other  hand,  it 
should  be  longer  than  the  camera  readout  time  to  allow  for  continuous 
imaging.  For  highly  radioactive  samples,  these  2  conditions  cannot  be 
simultaneously  satisfied;  therefore,  the  exposure  should  be  set  to  the 
highest  value  that  causes  no  overlap  between  tracks.  Typically,  this 
value  corresponds  to  an  average  of  10  ionization  events  per  frame 
(Supplemental  Fig.  2),  but  it  may  vary  with  the  spatial  distribution 
of  the  activity.  Loss  of  sensitivity  due  to  slow  camera  readout  is 
spatially  uniform.  Given  the  range  of  camera  exposures  available 
(10  ps-10  s),  the  microscope  can  image  a  wide  range  of  radionuclide 
activities,  from  0.2  Bq/mm2  to  20  kBq/mm2.  For  a  fixed  exposure 
time,  the  event  rate  is  also  linear  with  activity  concentration  (Supple¬ 
mental  Fig.  3). 

The  automated  extraction  of  (3  trajectories  from  the  sequence  of 
short  frames  involves  a  certain  number  of  image  processing  steps. 
Briefly,  each  frame  is  first  filtered  with  a  gaussian  kernel  to  reduce 
spatially  uncorrelated  shot  noise.  An  H-maxima  transform  is  then 
applied  to  suppress  small  local  image  maxima.  The  processed  image 
is  later  segmented  using  a  constant  threshold  set  above  the  noise  floor. 
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Morphologic  operations  are  also  applied  to  regularize  the  delineation 
of  the  segmented  tracks. 

For  image  processing  purposes,  we  distinguish  3  types  of  ionization 
tracks.  Long  tracks  (<500  pm)  are  produced  by  positrons  with  high 
initial  kinetic  energy  and  are  traced  backward  to  the  decay  origin 
using  a  custom  algorithm  (Fig.  2A).  Because  the  direction  of  the  par¬ 
ticle  along  the  track  is  ambiguous,  the  origin  is  defined  as  the  end  of 
the  track  closest  to  a  cell  (segmented  from  a  bright-field  micrograph). 
Shorter  tracks  are  characteristic  of  positrons  that  either  travel  a  short 
distance  because  of  lower  kinetic  energy  or  travel  orthogonally  to  the 
imaging  plane  (Fig.  2B).  These  tracks  are  directly  assigned  to  the 
nearest  cell.  Last,  annihilation  photons  are  more  likely  to  interact  deep 
in  the  scintillator  plate,  away  from  the  focal  plane  (Fig.  1A).  These 
interactions  sometimes  appear  as  a  diffuse,  out-of-focus  luminescence 
signal  and  are  rejected  during  image  processing  (Fig.  2C).  A  more 
detailed  description  of  the  track  reconstruction  algorithm  is  presented 
as  supplemental  information. 

The  final  ORBIT  image  is  reconstructed  by  aggregating  the  decay 
locations  estimated  from  the  sequence  of  frames.  The  resulting  image 
is  then  filtered  with  a  gaussian  kernel  to  account  for  the  localization 
uncertainty. 
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FIGURE  2.  Examples  of  ionization  tracks.  (A)  Long  ionization 
track,  originating  from  cell.  Dashed  line  is  estimated  positron  trajec¬ 
tory.  (B)  Shorter  ionization  track,  likely  produced  by  low-energy 
positron  or  positron  traveling  orthogonally  to  imaging  plane.  Red 
dot  is  estimated  decay  location.  (C)  Annihilation  photon  interacting 
deep  in  scintillator  and  producing  diffuse,  unfocused  signal.  Such 
signals  are  not  used  in  final  image.  Number  of  photoelectrons  (pe) 
recorded  by  CCD  sensor  for  each  event  is  indicated. 


Annihilation  and  Bremsstrahlung  Photon  Contamination 

To  evaluate  the  proportion  of  events  that  is  produced  by  secondary 
annihilation  and  Bremsstrahlung  photons,  we  imaged  the  ionization 
tracks  produced  by  a  thin  layer  of  18F-FDG  (19  kBq/mm2)  while  vary¬ 
ing  the  focal  plane  depth.  For  each  depth,  we  counted  how  many 
ionization  tracks  were  in  focus  in  200  sequential  images  (27-ms  expo¬ 
sure).  Because  the  objective  has  a  shallow  depth  of  field,  only  ioniza¬ 
tion  tracks  near  the  focal  plane  are  in  focus.  Because  of  the  limited 
range  of  positrons,  tracks  observed  deep  in  the  scintillator  are  attrib¬ 
uted  to  high-energy  photons. 

Cell  Imaging 

As  a  demonstration,  we  imaged  the  uptake  of  18F-FDG  in  human 
breast  adenocarcinoma  cells  and  compared  it  with  an  optical  fluores¬ 
cent  glucose  analog,  2-[A-(7-nitrobenz-2-oxa-l,3-diaxol-4-yl)amino]- 
2-deoxyglucose  (2-NBDG).  We  seeded  104  MDA-MB-231  cells  on  a 
10  x  10  x  0.5  mm  CdW04  scintillator  plate.  The  plate  was  immersed 
in  RPMI  1640  supplemented  with  10%  fetal  bovine  serum  and  1% 
penicillin/streptomycin  and  stored  at  37°C  in  a  humidified  atmosphere 
with  5%  C02.  The  next  day,  the  cells  were  approximately  10%  con¬ 
fluent  and  did  not  present  signs  of  altered  function  (Supplemental  Fig. 
4).  After  keeping  the  cells  fasting  for  30  min  in  glucose-free  RPMI, 
we  incubated  them  for  1  h  at  37°C  with  18F-FDG  (11  MBq)  and  2- 
NBDG  (100  pM)  and  washed  them  to  remove  the  excess  probe. 

Imaging  was  performed  by  2  methods.  First,  a  single  radiolumi¬ 
nescence  image  was  acquired  with  an  exposure  time  of  5  min  and  an 
electron-multiplication  gain  of  201.  Then,  a  sequence  of  6,000  frames 
was  acquired  with  an  exposure  time  of  50  ms — effectively  spanning 
5  min — and  an  electron-multiplication  gain  of  1,200  and  reconstructed 
into  a  single  image  using  ORBIT.  Image  reconstruction  took  about  8  h 
using  MATLAB  (The  Math  Works)  on  a  recent  computer,  with  most  of 
the  computation  spent  reconstructing  long  positron  tracks. 

RESULTS 

The  point- spread  function  and  modulation  transfer  function 
were  found  to  depend  strongly  on  the  depth  of  the  source  (Sup¬ 
plemental  Fig.  5).  The  Nyquist  frequency  for  the  camera  (0.15 
cycles/pm  for  2  x  2  binning)  was  attenuated  by  50%  at  a  2-pm 
depth  and  by  90%  at  4  pm.  For  a  small  point- source  object,  the 
response  of  the  microscope  decreased  gradually  with  depth  (Fig. 
3),  an  effect  known  as  optical  sectioning.  It  follows  from  these 
results  that  the  sensitive  region  of  the  scintillator  spans  approxi¬ 
mately  5  pm  on  either  side  of  the  focal  plane.  Therefore,  when  an 
energetic  positron  intersects  the  focal  plane,  the  microscope  can 
capture  only  the  portion  of  its  track  that  is  in  this  sensitive  region. 
The  simulation  of  simple  positron  trajectories  showed  that  the 
transition  between  the  2  regions  is  gradual  (Supplemental  Fig. 
6).  Because  of  optical  sectioning,  the  microscope  captures  only 
a  fraction  of  the  total  energy  deposited  by  energetic  positrons, 
which  can  travel  up  to  500  pm  in  CdW04. 

By  varying  the  depth  of  the  focal  plane,  we  estimated  that 
positrons  constitute  92%  of  the  events  recorded  at  the  scintillator 
edge  (Fig.  4).  The  event  counting  rate  was  1,010  cps  for  positrons 
and  90  cps  for  annihilation  and  Bremsstrahlung  photons.  Posi¬ 
trons  rarely  traveled  more  than  20  pm  from  the  scintillator  edge, 
highlighting  the  need  for  precise  focusing  to  maximize  positron 
sensitivity. 

An  ORBIT  radioluminescence  image  of  cellular  18F-FDG  up¬ 
take  was  synthesized  from  33,253  events  (110.8  cps).  As  a  control, 
we  also  imaged  a  sample  not  incubated  with  18F-FDG  and  mea¬ 
sured  a  background  rate  of  1 .2  cps,  likely  caused  by  cosmic  radi¬ 
ation  and  background  noise  in  the  camera.  On  the  basis  of  the 


Radioluminescence  Microscopy 


Pratx  et  al.  1843 


FIGURE  3.  Maximum  point-source  response  as  function  of  depth. 
Maximum  response  of  microscope  to  small  object  (gaussian  dis¬ 
tribution  with  a  full-width  half-maximum  of  6.5  |xm)  placed  at 
various  depths  was  computed  using  analytic  model.  Inset  shows 
2-dimensional  point-spread  function  (depth,  0,  1,  1.5,  and  2  |xm). 

camera  gain  and  calibration  factor,  we  estimated  that  single  ion¬ 
ization  events  created  between  20  and  1,000  photoelectrons  in  the 
CCD  sensor,  with  an  average  of  390  (Supplemental  Fig.  7A).  On 
the  basis  of  the  estimated  optical  sensitivity  of  the  system  and  the 
light  yield  of  the  scintillator,  these  correspond  to  minimum,  aver¬ 
age,  and  maximum  energy  depositions  of  10,  200,  and  510  keV, 
respectively.  These  observed  energy  depositions  reflect  the  diver¬ 
sity  of  positron  energies  and  trajectories  and  represent  only  ioni¬ 
zation  occurring  within  5  jim  of  the  focal  plane.  This  may  be  why 
the  average  measured  energy  was  smaller  than  the  average  posi¬ 
tron  energy  for  18F  (250  keV). 

For  both  imaging  modes,  the  noise  is  dominated  by  shot  noise 
and,  for  the  high-gain  acquisition,  excess  noise  factor,  with  a  neg¬ 
ligible  contribution  from  dark  noise  (0.001  electron/pixel/s)  and 
read  noise  (0.01  electron/pixel  at  11  MHz,  relative  to  electron- 
multiplication  gain).  On  average,  a  positron  produced  0.6  photo¬ 
electrons  per  CCD  pixel  (before  electron-multiplication  and 
analog  gain,  but  after  2x2  pixel  binning;  Supplemental  Fig.  7B). 
In  comparison,  a  pixel  receives  approximately  150  photoelectrons 
for  a  5-min  exposure.  As  a  result,  single  ORBIT  frames  have  a  low 
signal-to-noise  ratio  compared  with  longer  exposures  (Supplemen¬ 
tal  Table  1).  Other  sources  of  background  signal  include  Cerenkov 
luminescence — approximately  1.4  photons  per  decay  for  18F 
(13) — and  scintillator  afterglow  (e.g.,  lutetium  orthosilicate). 

Although  isolated  cells  are  clearly  distinguishable  on  the  5-min 
exposure  image  (Fig.  5A,  red  arrows),  ORBIT  provides  better 
resolution  of  single  cells  in  small  clusters  of  2-4  cells  (Fig.  3 A, 
yellow  arrows).  A  line  profile  through  the  images  confirms  that 
peaks  and  valleys  are  better  defined  on  the  ORBIT  image  and 
coincide  to  a  greater  degree  with  the  fluorescence  image  (Fig.  5B; 
Pearson  correlation  coefficient,  0.85  vs.  0.67).  In  contrast,  the  5- 
min-exposure  radioluminescence  image  presents  a  higher  level  of 
background  signal  between  peaks,  likely  caused  by  annihilation 
and  Bremsstrahlung  photons,  and  a  widening  of  the  peaks,  due  to 
positron  range  blurring.  In  addition,  ORBIT  directly  counts  how 
many  decays  occur  within  each  cell,  with  no  dependence  on  the 
variable  positron  energy.  Thus,  a  region-of-interest  analysis  found 
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that,  on  a  cell-by-cell  basis,  2-NBDG  fluorescence  was  more 
highly  correlated  with  18F-FDG  uptake  measured  with  ORBIT 
than  with  a  5-min  exposure  (Fig.  5C).  A  perfect  correlation  be¬ 
tween  18F-FDG  and  2-NBDG  is  not  achievable  as  the  2  probes  are 
not  strictly  identical:  2-NBDG  is  degraded  after  cell  internaliza¬ 
tion,  becoming  inactive,  whereas  18F-FDG  accumulates  in  the  cell 
as  FDG-6-phosphate.  Physical  differences  between  radionuclide 
decay  and  fluorescence  emission  can  also  account  for  discrepan¬ 
cies  between  the  2  signals. 

DISCUSSION 

We  found  that  ORBIT  could  visualize  individual  decay  events 
and  resulting  ionization  tracks.  By  processing  those  ionization 
tracks,  we  produced  images  of  the  distribution  of  18F-FDG  in 
living  cells  with  spatial  resolution  greater  than  the  positron  range. 
Whereas  a  single-exposure  radioluminescence  image  has  arbitrary 
camera  units,  an  ORBIT  image  has  units  of  counting  rate.  In  prin¬ 
ciple,  a  particle-counting  approach  such  as  ORBIT  should  be  more 
quantitative  than  an  integrating  approach.  That  may  be  why  OR¬ 
BIT  images  provided  measurements  of  18F-FDG  uptake  in  single 
cells  that  presented  a  stronger  correlation  to  the  uptake  of  a  fluo¬ 
rescent  glucose  analog.  A  precise  comparison  between  the  2  ap¬ 
proaches  will  be  the  aim  of  future  work. 

We  also  found  that  because  of  the  shallow  depth  of  field  of  the 
objective,  imaging  of  ionizing  events  is  restricted  to  a  10-|jim- thick 
slab  centered  on  the  focal  plane.  This  feature  of  the  microscope 
allowed  us  to  better  discriminate  between  positrons,  which  travel 
short  distances,  and  secondary  photons,  which  interact  more  deeply 
in  the  scintillator. 

ORBIT  works  best  for  nonconfluent  cell  cultures  because  each 
positron  track  can  reliably  be  assigned  to  a  cell  delineated  from 
a  bright-held  microscopy  image.  ORBIT  is  generally  able  to  dis¬ 
tinguish  individual  cells  in  small  clusters  (2-4  cells)  that  naturally 
occur  during  cell  growth  (Fig.  3A,  yellow  arrows).  However,  for  more 
confluent  cell  distributions  such  as  tissue  sections,  the  bright-held 
image  is  less  useful  for  estimating  the  positron  direction  along  the 
track  because  more  than  one  cell  may  be  equally  close  to  the 


FIGURE  4.  Annihilation  and  Bremsstrahlung  photon  contamina¬ 
tion.  A  thin  layer  of  18F-FDG  (19  kBq/mm2)  placed  between  scintil¬ 
lator  and  imaging  dish  was  imaged  by  acquiring  200  frames  (3-ms 
exposure).  Events  in  focus — that  is,  that  displayed  high-frequency 
features — were  counted  at  various  focal  depths. 
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FIGURE  5.  Radioluminescence  microscopy  of  18F-FDG  uptake.  Human  breast  cancer 
cells  (MDA-MB-231)  were  deprived  of  glucose  for  1  h,  incubated  for  1  h  with  18F-FDG 
(11  MBq)  and  2-NBDG  (100  piM),  and  then  washed.  (A)  Bright-field  (scale  bar,  200  |xm), 
fluorescence  (2-NBDG),  single-exposure  radioluminescence  microscopy  (FDG-RM),  and 
ORBIT  (FDG-ORBIT).  Objective  is  40x/1.3  numerical  aperture.  (B)  Radioluminescence  and 
fluorescence  intensity  shown  along  line  profile  (red  dashed  line  in  A).  (C)  Correlation  be¬ 
tween  cellular  uptake  of  2-NBDG  and  18F-FDG,  where  18F-FDG  uptake  was  measured  with 
conventional  radioluminescence  microscopy  (FDG-RM)  or  ORBIT  (FDG-ORBIT).  Cellular 
uptake  was  computed  by  drawing  circular  region  of  interest  around  each  cell  on  bright-field 
image. 


ionization  track.  This  is  unfortunate  because  dense  cell  distribu¬ 
tions  are  more  representative  of  physiologic  conditions.  ORBIT 
can  be  adapted  to  imaging  dense  tissue  in  several  ways.  As  the 
positron  slows  down,  it  becomes  more  ionizing  and  its  ionization 
track  becomes  brighter.  Furthermore,  as  the  positron  moves  away 
from  the  focal  plane,  the  ionization  track  becomes  defocused. 
These  properties  may  help  determine  the  direction  of  the  positron 
along  its  track.  Alternatively,  long  ionization  tracks  may  also  be 
discarded  in  favor  of  shorter  tracks  that  are  easier  to  localize. 
Tissue  dissociation  techniques  may  also  allow  dense  tissues  to 
be  physically  separated  for  imaging  (14). 

The  ORBIT  method  can  readily  be  extended  to  dynamic  studies. 
Short  frames  can  be  acquired  for  the  desired  duration  of  the  study, 
retrospectively  gated  into  time  intervals,  and  independently  re¬ 
constructed.  However,  the  use  of  the  bright-held  image  to  de¬ 
lineate  cells  is  not  optimal  in  the  presence  of  background  activity. 
Positrons  emitted  from  a  free  radionuclide  probe  will  produce 
ionization  tracks  that  should  not  be  associated  with  a  cell.  Thus, 
to  avoid  bias,  p  particles  should  be  localized  without  using 
the  bright-held  image  when  signihcant  activity  is  present  in  the 
background. 

Highly  sensitive  CCD  cameras  with  high  internal  gain,  com¬ 
parable  to  the  one  used  in  this  work,  have  in  recent  years  found 
application  in  high-resolution  SPECT  imaging  (15-18).  These 
SPECT  systems  image  the  ionization  track  produced  by  low-en¬ 
ergy  7  rays  interacting  in  a  thin  scintillator  screen.  Similar  to  our 
work,  they  apply  image-processing  algorithms  to  extract  the  po¬ 
sition  and  energy  of  each  individual  y-ray  interaction.  Such  an 
imaging  scheme  enables  spatial  resolution  on  the  order  of  50-70 
|jim,  much  higher  than  conventional  SPECT  systems.  In  a  similar 
fashion,  a  digital  autoradiography  system  was  recently  developed 
using  a  high-efficiency  CCD  sensor  (19).  By  placing  the  tissue  of 


interest  in  direct  contact  with  the  CCD, 
this  system  can  directly  visualize  the  tra¬ 
jectories  of  positrons  entering  the  silicon 
detector.  The  system  achieved  a  spatial 
resolution  of  86  jim  when  imaging  18F  in 
a  100- pan-thick  tissue  section.  However,  it 
is  not  compatible  with  live  cell  imaging 
because  the  CCD  sensor  must  be  cooled 
to  very  low  temperatures.  Furthermore,  un¬ 
like  radioluminescence  microscopy,  it  is 
unable  to  image  the  tissue  optically. 


CONCLUSION 

Because  of  limited  spatial  resolution, 
autoradiography  and  other  methods  can 
measure  the  uptake  of  a  radiotracer  only  in 
a  relatively  large  cell  population.  In  con¬ 
trast,  radioluminescence  microscopy  is 
able  to  measure  radiotracer  uptake  with 
single-cell  resolution,  in  the  familiar  en¬ 
vironment  of  a  fluorescence  microscope, 
using  fluorescent  probes  to  highlight  bio¬ 
chemical  features  of  interest.  Dynamic 
studies  can  also  shed  light  on  radiotracer 
kinetics  in  vitro  and  at  the  single-cell  level. 
The  ORBIT  scheme  presented  here  makes 
those  single-cell  measurements  more  ac¬ 
curate  with  no  increase  in  imaging  time  or 
hardware  complexity.  Future  work  will  focus  on  applying  this 
approach  for  imaging  denser  cell  distributions  such  as  tissue 
sections. 
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Erratum 

In  the  article  “Assessment  of  Treatment  Response  in  Patients  with  Glioblastoma  Using  G-(2-18F-Fluoroethyl)-L- 
Tyrosine  PET  in  Comparison  to  MRI,”  by  Galldiks  et  al.  (J  Nucl  Med.  2012;53:1048-1057),  the  abstract  mistakenly 
stated,  “Receiver-operating-characteristic  analysis  showed  that  a  decrease  of  the  TBR(max)  between  FET-1  and 

FET-2  of  more  than  20%  predicted  poor  survival - ”  However,  such  a  decrease  in  fact  predicted  favorable  survival, 

not  poor  survival.  The  authors  regret  the  error. 
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Cerenkov  luminescence  imaging  (CLI)  is  an  emerging  new 
molecular  imaging  modality  that  is  relatively  inexpensive, 
easy  to  use,  and  has  high  throughput.  CLI  can  image 
clinically  available  PET  and  SPECT  probes  using  optical 
instrumentation.  Cerenkov  luminescence  endoscopy  (CLE)  is 
one  of  the  most  intriguing  applications  that  promise  potential 
clinical  translation.  We  developed  a  prototype  customized 
fiberscopic  Cerenkov  imaging  system  to  investigate  the 
potential  in  guiding  minimally  invasive  surgical  resection. 
Methods:  All  experiments  were  performed  in  a  dark  cham¬ 
ber.  Cerenkov  luminescence  from  18F-FDG  samples  contain¬ 
ing  decaying  radioactivity  was  transmitted  through  an  optical 
fiber  bundle  and  imaged  by  an  intensified  charge-coupled 
device  camera.  Phantoms  filled  with  18F-FDG  were  used  to 
assess  the  imaging  spatial  resolution.  Finally,  mice  bearing 
subcutaneous  C6  glioma  cells  were  injected  intravenously 
with  18F-FDG  to  determine  the  feasibility  of  in  vivo  imaging. 
The  tumor  tissues  were  exposed,  and  CLI  was  performed  on 
the  mouse  before  and  after  surgical  removal  of  the  tumor 
using  the  fiber-based  imaging  system  and  compared  with 
a  commercial  optical  imaging  system.  Results:  The  sensitiv¬ 
ity  of  this  particular  setup  was  approximately  45  kBq  (1.21 
|jiCi)/300  |ulL.  The  3  smallest  sets  of  cylindric  holes  in  a  com¬ 
mercial  SPECT  phantom  were  identifiable  via  this  system, 
demonstrating  that  the  system  has  a  resolution  better  than 
1.2  mm.  Finally,  the  in  vivo  tumor  imaging  study  demon¬ 
strated  the  feasibility  of  using  CLI  to  guide  the  resection  of 
tumor  tissues.  Conclusion:  This  proof-of-concept  study  ex¬ 
plored  the  feasibility  of  using  fiber-based  CLE  for  the  detection 
of  tumor  tissue  in  vivo  for  guided  surgery.  With  further  improve¬ 
ments  of  the  imaging  sensitivity  and  spatial  resolution  of  the  cur¬ 
rent  system,  CLE  may  have  a  significant  application  in  the  clinical 
setting  in  the  near  future. 
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V^erenkov  luminescence  imaging  (CLI)  has  recently 
attracted  increasing  interest  in  the  field  of  molecular  im¬ 
aging  (7,2).  CLI  is  a  new  optical  imaging  modality  in 
which  images  are  obtained  by  monitoring  the  Cerenkov 
photons  emitted  from  highly  energetic  moving  charged 
particles  (|3+  or  |3_).  Conventional  nuclear  imaging 
methods,  such  as  PET  and  SPECT,  are  the  most  widely 
used  clinical  molecular  imaging  techniques.  However, 
these  modalities  usually  suffer  from  high  cost,  limited 
availability,  relatively  low  spatial  resolution,  and  low 
throughput  (J).  As  a  technique  that  bridges  optical  imaging 
and  radionuclide  imaging,  CLI  has  shown  many  advantages 
such  as  high  sensitivity,  high  resolution,  low  cost,  wide  avail¬ 
ability,  relatively  high  throughput,  and  commercially  avail¬ 
able  radionuclide  probes  already  approved  by  the  Food  and 
Drug  Administration. 

Since  its  discovery  in  2009,  CLI  has  quickly  become 
a  practical  molecular  imaging  technique,  and  many  new 
applications  of  CLI  in  preclinical  research  continue  to 
emerge  (7,2).  Several  research  groups  have  demonstrated 
that  CLI  can  be  a  powerful  tool  for  tumor  imaging  using 
radionuclide  probes  such  as  18F-FDG.  Important  validation 
studies  have  also  been  performed,  and  it  was  reported  that 
there  is  a  good  linear  correlation  between  the  tumor  uptake 
quantified  by  PET  and  tumor  CLI  signals  in  subcutaneous 
xenograft  models  (4-7).  CLI  has  also  been  found  to 
be  useful  for  imaging  a  and  pure  |3“  emitters  such  as  90 Y 
and  225 Ac,  which  are  used  for  cancer  treatment  (4).  More¬ 
over,  CLI  can  be  applied  in  the  monitoring  of  reporter  gene 
expression.  The  herpes  simplex  virus  type  1  thymidine  ki¬ 
nase  and  the  sodium  iodide  symporter  reporter  genes  were 
recently  successfully  imaged  by  CLI  in  conjunction  with 
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appropriate  radioactive  reporter  probes  {4,8,9).  Cerenkov 
photons  have  also  been  used  to  illuminate  fluorescent 
dyes  and  nanoparticles  for  in  vivo  imaging  (10-13). 
A  recent  study  further  showed  that  CLI  could  be  used 
for  monitoring  tumor  surgery  in  xenograft  mouse  models, 
highlighting  the  translational  potential  of  the  modality 
(14).  CLI  has  also  been  used  to  measure  the  radiochem¬ 
ical  purity  of  a  radiolabeled  compound  and  to  image 
plant  physiology  (75).  In  addition  to  these  applications, 
the  recent  advancement  in  CLI  tomography  makes  this 
novel  imaging  modality  even  more  powerful  and  prom¬ 
ising  for  biomedical  research  including  diagnostic  imag¬ 
ing  and  therapeutic  monitoring  (16-19). 

Because  it  can  image  clinically  available  radiotracers, 
CLI  has  the  potential  to  be  rapidly  translated  into  clinical 
applications  (4,18).  However,  all  the  studies  described 
above  have  involved  the  use  of  conventional  small-animal 
optical  imaging  systems,  which  are  not  compatible  with 
routine  clinical  practice.  In  this  study,  we  have  built  the 
first,  to  our  knowledge,  prototype  system  that  is  amenable 
to  Cerenkov  luminescence  endoscopy  (CLE)  in  the 
clinic.  This  system  comprises  an  optical  fiber  bundle 
and  an  intensified  charge-coupled  device  (CCD)  camera. 
The  small-diameter  flexible  endoscope  is  designed  for  min¬ 
imally  invasive  monitoring  of  living  tissues  and  organs, 
using  a  sensitive  camera  to  produce  images  of  the  radionu¬ 
clide  probe.  It  is  also  expected  that  hollow  organs  (e.g., 
bladder  and  lung)  or  insufflated  cavities  (e.g.,  esophagus 
and  colon)  inside  the  body  can  provide  natural  and  ana¬ 
tomic  dark  chambers  for  CLE.  By  overlaying  conventional 
bright-field  images  commonly  obtained  from  an  endoscope 
with  the  Cerenkov  luminescence  images  generated  by  ra¬ 
dionuclide  probes,  CLE  could  be  used  to  identify  diseased 
tissues  for  diagnostic  purposes  and  real-time  monitoring  of 
endoscopic  surgery. 

Therefore,  we  investigated  the  feasibility  of  CLE  for 
guiding  cancer  surgical  resection  by  performing  sensi¬ 
tivity,  spatial-resolution,  and  proof-of-concept  experi¬ 
ments.  The  characteristic  sensitivity  of  this  system  was 
evaluated  using  the  natural  decay  of  18F-FDG.  Imaging 
resolution  was  determined  using  a  phantom  with  cylindric 
holes.  Finally,  imaging  studies  were  performed  to  dem¬ 
onstrate  surgical  guidance  in  a  small-animal  tumor  mo¬ 
del.  A  well-known  commercial  optical  imaging  system 
was  used  to  provide  a  comparison.  The  data  presented 
herein  outline  the  instrumentation  and  methodology;  fur¬ 
ther  steps  that  are  needed  to  translate  this  modality  into 
the  clinic  are  also  discussed. 

MATERIALS  AND  METHODS 

18F-FDG  was  produced  by  the  Radiochemistry  Facility  at  Stan¬ 
ford  University.  The  rat  glioma  cell  line  C6  was  obtained  from 
American  Type  Culture  Collection.  Female  athymic  nude  mice 
(nu/nu)  (n  =  5)  obtained  from  Charles  River  Laboratories,  Inc., 
were  4-6  wk  old.  A  CRC-15R  PET  dose  calibrator  (Capintec  Inc.) 
was  used  for  all  radioactivity  measurements. 


Tumor  Model 

All  animal  studies  were  performed  in  compliance  with  federal 
and  local  institutional  guidelines  for  the  conduct  of  animal 
experimentation.  C6  cells  were  cultured  in  Dulbecco  modified 
Eagle  medium  supplemented  with  10%  fetal  bovine  serum  and  1% 
penicillin/streptomycin  (Invitrogen  Life  Technologies).  The  cell 
line  was  maintained  in  a  humidified  atmosphere  of  5%  C02  at 
37°C,  with  the  medium  changed  every  other  day.  A  75%  confluent 
monolayer  was  detached  with  trypsin  and  dissociated  into  a  sin¬ 
gle-cell  suspension  for  further  cell  culture.  Approximately  1  x  106 
C6  cells  suspended  in  phosphate-buffered  saline  (0.1  M,  pH  7.2; 
Invitrogen)  were  implanted  subcutaneously  in  the  left  legs  of  nude 
mice.  Tumors  were  allowed  to  grow  to  150-200  mm3  (2-3  wk), 
and  the  tumor-bearing  mice  were  imaged  in  vivo. 

CLI 

Validation  CLI  was  performed  with  an  IVIS  Spectmm  system 
(Caliper  Life  Sciences).  For  all  in  vivo  studies,  radionuclides  were 
diluted  in  phosphate-buffered  saline.  Animals  were  placed  in 
a  light-tight  chamber  under  isoflurane  anesthesia.  Each  acquisition 
took  3  min  for  all  studies  without  filters.  Images  were  acquired  and 
analyzed  using  Living  Image  3.0  software  (Caliper  Life  Sciences). 
The  mice  were  kept  fasting  overnight  before  18F-FDG  imaging  and 
anesthetized  by  inhalation  of  2%  isoflurane  during  the  experiment. 

Fiber-Based  CLE  System 

The  fiber-based  CLE  system  used  an  imaging  optical  fiber 
coupled  to  a  highly  sensitive  intensified  CCD  camera.  The 
application  of  this  system  for  surgical  oncology  guidance  is 
detailed  in  Figure  1.  Specifically,  a  microimaging  lens  (Cinegon, 
F/1.4,  12-mm  focal  length;  Schneider)  was  coupled  at  the  distal 
end  of  an  optical  imaging  fiber  bundle  that  was  108  mm  long,  with 
a  5  x  6.7  mm  active  area.  The  pixels  of  the  fiber  bundle  are  made 
of  discrete  10- pm  fibers  (IG154;  Schott).  At  the  proximal  end, 
a  Pentax  (F/1.4)  lens  provided  relay  optics  to  the  camera  (Supple¬ 
mental  Fig.  1A;  supplemental  materials  are  available  online  only 
at  http://jnm.snmjoumals.org).  The  camera,  an  image-intensified 


FIGURE  1.  Suggested  application  of  fiber-based  system  for  en¬ 
doscopic  and  laparoscopic  CLI. 
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CCD  (Turbo  640-Z;  Stanford  Photonics  Inc.),  had  640  x  480 
pixels  and  single-photon  imaging  capability  (Supplemental  Fig. 
IB).  To  minimize  background  light  and  to  emulate  a  light-tight 
anatomic  cavity,  all  images  were  taken  in  a  dark  box.  In  addition, 
images  were  postprocessed  with  a  software  thresholding  tech¬ 
nique  that  removed  pixel  values  below  a  manufacturer-recom¬ 
mended  value;  this  optimal  value  was  chosen  to  remove  low- 
intensity  pixels  that  corresponded  to  thermal  and  read  noise  on 
the  CCD.  These  steps  reduced  the  noise  to  approximately  20 
counts  per  second. 

Imaging  Using  Fiber-Based  System 

Images  were  acquired  at  a  high  frame -rate  of  120  Hz  to  allow 
for  maximal  reduction  in  noise.  The  noise-reduction  step  de¬ 
scribed  above  was  performed  for  each  frame.  These  frames  were 
then  accumulated  for  5  min  to  form  the  raw  data.  The  raw  data 
were  image-processed  offline  using  an  algorithm  to  remove 
cosmic  and  stray  y-events.  This  algorithm  removed  noisy  pixels 
by  thresholding  the  gradient  around  each  pixel  in  its  local  area.  An 
additional  denoising  step  removed  spurious  bright  pixels  that 
appeared  in  sequential  exposures. 

The  fiber-based  CLE  system  was  characterized  for  sensitivity 
by  imaging  2  adjacent  wells  of  a  black  96- well  plate  (300  pL  per 
well).  A  glycerol  and  water  mixture  combined  with  3.7  MBq  (100 
pCi)  of  18F-FDG  was  in  1  well,  and  a  glycerol  and  water  mixture 
only  was  in  a  nearby  well.  Images  were  sequentially  acquired  over 
approximately  10  half-lives  (19  h,  43  min),  and  each  data  point 
was  formed  from  5  min  of  photon  accumulation. 

The  spatial  resolution  of  the  system  was  characterized  using 
a  standard  PET/SPECT  phantom  (Micro  Hot-Spot  Phantom; 
Data  Spectrum  Corp.)  filled  with  15.2  MBq  (410  pCi)  of  18F- 
FDG;  this  phantom  is  also  commonly  used  to  evaluate  the  spatial 
resolution  of  PET  systems.  Both  the  2.4-  and  the  1.6-mm  cylin- 
dric  holes  were  used  to  form  line  profiles,  taken  from  the  fiber- 
based  images.  An  ambient-light  image  and  a  functional  Ceren¬ 
kov  luminescence  image  were  acquired.  The  exposure  time  for 
all  ambient  images  was  less  than  1  s,  and  the  functional  image 
acquisition  time  was  5  min.  The  subject  was  5  cm  from  the 
optical  system. 

Fiber-Based  CLE  System  for  Surgery  Monitoring 

Five  mice  were  injected  with  37  MBq  (1  mCi)  of  18F-FDG  via 
the  tail  vein.  The  tracer  was  allowed  to  accumulate  for  60-70  min, 
and  the  mice  were  imaged  in  the  I  VIS  system  to  verify  tumor 
uptake  of  18F-FDG.  Tumors  were  then  imaged  with  the  IVIS 
and  fiber-based  CLE  systems  after  removal  of  the  skin  covering 
the  subcutaneous  tumor  and  after  excision  of  the  tumor.  To  sim¬ 
ulate  an  environment  mimicking  surgical  resection,  and  to  validate 
the  origin  of  the  signal,  the  excised  tumor  was  also  imaged  after 
being  placed  adjacent  to  the  surgical  cavity.  Two  images  were 
acquired  for  each  of  these  steps:  1  ambient  image  and  1  functional 
Cerenkov  luminescence  image.  The  exposure  time  for  all  ambient 
images  was  less  than  1  s,  and  the  functional  image  acquisition 
time  was  5  min. 

Statistical  Methods 

Quantitative  data  were  expressed  as  mean  ±  SD.  Means  were 
compared  using  the  Student  t  test.  A  95%  confidence  level  was 
chosen  to  determine  the  significance  between  groups,  with  P  val¬ 
ues  of  less  than  0.05  indicating  significant  differences. 


FIGURE  2.  System  sensitivity  via  sequential  imaging  of  3.7  MBq 
(100  pCi)  of  18F-FDG  for  approximately  20  h.  Optical  signals  from 
18F-FDG  and  control  sample  were  plotted. 

RESULTS 

System  Characterization 

Figure  2  shows  the  reduction  in  signal  during  18F-FDG 
decay  over  time  for  2  wells:  1  filled  with  18F-FDG  and  1 
filled  with  a  water  and  glycerol  mixture  only.  With  a  5 -min 
integration  time,  a  minimum  of  approximately  45  kBq 
(1.21  pCi)  of  activity  can  be  identified  as  different  from 
the  control  well  containing  the  water  and  glycerol  solution 
(signal-to-noise  ratio  >1).  There  was  a  decrease  in  signal 
in  the  control  well  due  to  stray  7-photons  from  the  18F-FDG 
well  interacting  with  the  optical  fiber  and  creating  scintil¬ 
lation  background  light  in  the  image. 

Visually,  the  photographic  and  Cerenkov  images  taken 
using  the  standard  IVIS  imaging  system  (Supplemental  Fig. 
2)  are  similar  to  the  photographic  and  Cerenkov  images 
taken  with  the  fiber-based  system  (Fig.  3A).  The  line  pro¬ 
files  demonstrate  a  high  correlation  between  the  Cerenkov 
and  photographic  line  profiles.  For  the  2.4-mm  holes  (Fig. 
3B),  the  line  profile  limits  depicted  by  the  large  green  stars 
in  Figure  3 A  had  a  Pearson  correlation  coefficient  of  0.78 
(significance,  P  <  le— 5).  For  the  1.6-mm  holes  (Supple¬ 
mental  Fig.  3),  the  line  profile  limits  depicted  by  the  smaller 
red  stars  had  a  Pearson  correlation  coefficient  of  0.71  (sig¬ 
nificance,  P  <  le— 5).  The  peaks  were  also  in  near-identical 
locations;  the  difference  in  peak  locations  between  the  holes 
for  the  line  profiles  of  the  photographic  and  Cerenkov-emis- 
sion  images  was  3.8%,  a  difference  of  only  a  few  pixels. 

System  Demonstration 

The  images  comparing  both  systems  before  removal  of 
the  tumor  for  mouse  1  are  shown  in  Figure  4  (IVIS  system 
[Fig.  4 A]  and  fiber-based  CLE  system  [Fig.  4B]).  Figure  5 A 
shows  the  images  produced  by  the  IVIS  system,  and  Figure 
5B  shows  images  from  the  fiber-based  CLE  system  after 
removal  of  the  tumor  in  mouse  1. 

Comparing  images  in  Figure  5B  quantitatively,  we  de¬ 
termined  the  tumor-to-background  ratio  for  the  Cerenkov 
signal  by  computing  the  ratio  of  the  median  value  in  the 
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FIGURE  3.  Characterization  of  fiberscopic 
system  spatial  resolution.  (A)  Respective 
ambient  (left)  and  Cerenkov  images  (right) 
of  PET/SPECT  phantom  (4.3-cm  inner  diam¬ 
eter).  (B)  Quantitative  line  profiles  of  ambient 
and  Cerenkov  images  for  2.4-mm  holes;  line 
sampled  is  indicated  by  green  stars  in  A. 
a.u.  =  arbitrary  units. 
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region  of  interest  encircling  the  excised  tumor  to  the  me¬ 
dian  value  in  the  cleared  tumor  cavity.  This  ratio  was  1.28 
(for  comparison,  the  tumor-to-background  ratio  from  the 
IVIS  system  was  1.16  for  the  3-min  scan).  For  mice  2,  3, 
4,  and  5,  the  tumor-to-background  ratios  for  the  removed 
tumor  were  1.41,  1.21,  1.02,  and  1.17,  respectively.  Tumor 
tissue  light  emission  was  significantly  higher  than  the  exposed 
cavity  for  all  mice  (Student  t  test  for  paired  samples,  P  <  0.05 
for  all).  Much  of  the  low  tumor-to-background  value  in  mouse 
4  can  be  explained  by  the  reflection  of  the  light  emitted  from 
the  tumor  by  the  tumor  cavity,  which  was  directly  adjacent. 
Residual  tumor  tissue  after  surgery  is  also  possible. 

DISCUSSION 

In  this  study,  we  demonstrated  the  methodology  and 
performance  of  an  optical  fiber  system  built  to  image  tumor 
margins  during  surgical  resection.  Cerenkov  imaging  has 
been  proposed  for  small-animal  drug  discovery  and  surgical 
guidance,  yet  previous  studies  have  used  an  imaging  box 
with  a  lens-mounted  imaging  system;  this  setup  is  not  realistic 
for  subjects  that  are  larger  (e.g.,  large  animals  and  humans). 
We  recently  demonstrated  the  feasibility  of  endoscopic 
imaging  of  Cerenkov  light  using  conventional  optical  fiber 
bundle/clinical  endoscopes,  an  optical  imaging  lens  system, 
and  a  sensitive  low-noise  CCD  camera  (20).  Hereby,  for  the 
first  time  we  investigated  a  Cerenkov  endoscope  for  surgi¬ 
cal  resection.  We  have  shown  that  this  system  is  capable  of 
imaging  1.2-mm  structures,  using  phantom  experiments. 
We  determined  the  sensitivity  of  the  system  for  18F-FDG 
(—45  kBq  [1.21  [xCi]/300  p,L).  We  then  demonstrated  the 
ability  to  visualize  the  accumulation  of  18F-FDG  in  a  tumor 
using  this  system  and  could  track  the  contrast  as  the  tumor 
was  excised  from  the  mouse. 

A  benefit  of  a  Cerenkov  endoscope  is  that  it  provides 
higher  image  spatial  resolution  than  PET  and  SPECT 
cameras,  because  the  optical  photons  emitted  from  the 
tissue  surface  are  detected  with  optical  lenses  and  a  high- 
resolution  CCD  chip.  Cho  et  al.  demonstrated  a  limit  of 
350  |mm  (in  full  width  at  half  maximum)  with  a  system 
designed  for  visualizing  a  microfluidic  chip  setup  (27). 
This  resolution  is  affected  by  the  imaging  system,  the 
range  of  the  (3-particles  in  the  tissue,  and  tissue  optical 


scatter.  Levin  and  Hoffman  showed  that  the  mean  |3 -particle 
track  of  18F-FDG  had  a  full  width  at  half  maximum  value 
of  102  |mm  in  water  (22),  which  determines  the  physical 
resolution  limit  of  a  Cerenkov  scope  when  imaging  18F- 
FDG.  Improved  resolution  may  be  achieved  using  a  radio¬ 
isotope  with  a  (3-particle  decay  of  lower  energy,  such  as 
13 1I,  although  this  would  result  in  lower  sensitivity  be¬ 
cause  there  is  an  inverse  relationship  between  resolution 
and  sensitivity.  We  do  note  some  inconsistency  in  the 
image  of  the  1.2-mm  holes,  but  many  of  these  holes  are 
clearly  visualized.  The  high  spatial  resolution  of  the  sys¬ 
tem  requires  accurate  focus  of  the  optics  both  at  the  distal 
fiber-lens  junction  and  at  the  proximal  fiber-camera  junc¬ 
tion — an  optical  focus  that  will  be  improved  in  our  further 


FIGURE  4.  Mouse  1  bearing  C6  glioma  after  tail-vein  administra¬ 
tion  of  37  MBq  (1  mCi)  of  18F-FDG.  (A)  Mouse  was  imaged  by 
commercially  available  optical  IVIS  system,  and  images  were  com¬ 
pared  with  those  from  prototype  fiber-based  system  (B).  Tumor 
tissues  are  outlined  by  red  lines.  Ambient-light  images  are  on  left, 
luminescent  images  are  in  middle,  and  fused  images  are  on  right. 
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FIGURE  5.  Mouse  1  was  imaged  by  IVIS  optical  system  (A)  and 
fiber-based  system  (B)  after  surgery  to  remove  tumor  tissues. 
Ambient-light  images  are  on  left,  luminescent  images  are  in  middle, 
and  fused  images  are  on  right. 

studies.  In  the  case  of  our  phantom,  the  resolution  of  the 
Cerenkov  luminescence  images  was  also  degraded  by 
reflections  within  the  translucent  phantom  and  the  sig- 
nal-to-noise  ratio  of  the  single-photon-counting  camera, 
which  was  not  entirely  optimized.  The  signal-to-noise  ra¬ 
tio  of  our  ambient-light  image  was  also  affected  by 
a  slightly  imperfect  focus  and  the  suboptimal  settings  of 
the  illumination  source  and  camera.  Thus,  although  we 
demonstrated  the  ability  to  visualize  many  of  these  1.2- 
mm  cylindric  holes  in  a  phantom  with  this  endoscope- 
based  system,  we  expect  to  be  able  to  image  submillimeter 
lesions  with  a  more  optimized  system. 

The  sensitivity  limit  of  this  system  was  45  kBq  (1.2 
p,Ci)/300  p,L  with  18F-FDG.  To  put  this  in  a  practical  per¬ 
spective,  if  this  system  were  used  for  head  and  neck  tumors 
of  the  oropharynx  (with  a  standardized  uptake  value  of 
5.17),  2.1  GBq  (54  mCi)  would  have  to  be  injected  intra¬ 
venously  into  a  70-kg  patient  to  have  enough  activity  to 
visualize  a  300- pug  tumor  (23)  in  a  5 -min  scan;  a  more 
sensitive  scope  would  enable  the  identification  of  smaller 
structures  or  similar  structures  with  a  reduction  in  ionizing 
dose  to  the  patient  and  clinical  staff.  Increased  sensitivity 
may  be  realized  with  more  sensitive  optics  (such  as  an 
F/0.95  lens)  or  fiber  optic  glass,  which  transmits  farther 
into  the  violet  and  ultraviolet  emissions.  The  fiberglass  used 
in  this  experiment  was  Schott-75  glass,  which  transmits 
40%  light  at  500  nm;  thus,  it  is  not  optimized  for  detecting 
the  short-wavelength-dominant  Cerenkov  emission.  In  ad¬ 
dition,  7-photons  emitted  from  the  radiotracer  scintillate 
impurities  in  the  glass,  resulting  in  undesired  background 
noise.  The  use  of  fused  silica,  which  is  much  more  sensitive 


to  the  ultraviolet  and  violet  emissions  and  has  fewer  impu¬ 
rities,  would  improve  sensitivity. 

As  a  prototype  device,  the  off-the-shelf  optical  lens 
on  this  system  was  too  large  (3-cm  diameter)  for  many 
endoscopic  or  laparoscopic  applications.  However,  smaller 
custom  optics  could  be  designed  to  enable  imaging  through 
smaller  apertures.  Meanwhile,  optimized  optic  lenses  and 
fibers  can  provide  better  sensitivity  and  shorter  integra¬ 
tion  time  in  upcoming  systems.  We  leave  this  investiga¬ 
tion  for  future  studies  using  more  customized  optical 
systems. 

The  main  advantage  of  this  system  over  other  optical 
techniques,  such  as  fluorescence  imaging,  is  that  it  is  able 
to  use  standard  PET  tracers,  such  as  18F-FDG,  which  are 
already  available  in  the  clinic.  Many  novel  tracers  such  as 
3'-deoxy-3'-18F-fhiorothymidine,  engineered  proteins,  and 
labeled  arginine-glycine-aspartic  acid  peptides  are  in  the 
pipeline  of  clinical  trials  required  by  the  Food  and  Drug 
Administration  because  of  their  utility  in  PET.  This  is  a  sig¬ 
nificant  advantage  for  this  technique,  because  much  effort 
and  innovation  is  being  dedicated  to  PET  tracers.  In  con¬ 
trast,  because  of  the  lack  of  a  standardized  imaging  device 
for  optical  imaging,  optical  contrast  agents  are  not  being 
sufficiently  developed  to  keep  pace  with  radiotracers.  Thus, 
the  main  advantage  of  a  Cerenkov  emission  surgical  scope 
is  that  it  is  readily  clinically  translatable.  The  combination 
of  PET  and  Cerenkov  imaging  provides  an  opportunity  to 
use  an  identical  signal  origin,  and  perhaps  an  identical  in¬ 
jection,  to  localize  tumors  for  excision. 

CONCLUSION 

This  proof-of-concept  study  demonstrated  the  feasibility 
of  using  fiber-based  CLE  for  the  detection  of  tumor  tissue 
in  vivo  and  demonstrated  its  potential  use  for  image-guided 
surgery.  With  further  improvement  in  imaging  sensitivity  and 
resolution  of  the  current  system,  it  is  expected  that  Cerenkov 
imaging  might  soon  be  translated  into  clinical  applications. 
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